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Does the observable spectrum of cosmic fluctuations depend on detailed initial conditions? This 
addresses the question if the general inflationary paradigm is sufficient to predict within a given 
model the spectrum and amplitude of cosmic fluctuations, or if additional particular assumptions 
about the initial conditions are needed. The answer depends on the number of e-foldings Ni„ 
between the beginning of inflation and horizon crossing of the observable fluctuations. We discuss 
an interacting inflaton field in an arbitrary homogeneous and isotropic geometry, employing the 
quantum effective action T. An exact time evolution equation for the correlation function involves 
the second functional derivative The operator formalism and quantum vacua for interacting 

fields are not needed. Use of the effective action also allows one to address the change of frames 
by field transformations (field relativity). Within the approximation of a derivative expansion for 
the effective action we find the most general solution for the correlation function, including mixed 
quantum states. For not too large Wn the memory of the initial conditions is preserved. In this 
case the cosmic microwave background cannot disentangle between the initial spectrum and its 
processing at horizon crossing. The inflaton potential cannot be reconstructed without assumptions 
about the initial state of the universe. We argue that for very large Wn a universal scaling form of the 
correlation function is reached for the range of observable modes. This can be due to symmetrization 
and equilibration effects, not yet contained in our approximation, which drive the short distance 
tail of the correlation function towards the Lorentz invariant propagator in flat space. 
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Are amplitude and spectrum of the primordial cosmic 
fluctuations uniquely determined by the inflaton potential 
or Is there some Influence of the Initial state of the uni¬ 
verse? Is there a good reason to assume a quantum vacuum 
state during inflation and is this unique? Is it necessary to 
postulate or understand particular properties of the initial 
conditions in order to explain the simple successful picture 
of the almost scale invariant spectrum? Or could It be 
that the dynamics Is such that all memory of the correla¬ 
tion function about the initial state is lost, such that the 
two point function during inflation is essentially uniquely 
fixed for a given solution of the field equations? Answers to 
these questions are crucial in order to assess if the inflation¬ 
ary paradigm by itself is sufficient for a unique prediction 
of the properties of primordial fluctuations within a given 
model, or If additional ingredients are needed for their un¬ 
derstanding. 

Conceptually, It is rather obvious that the fluctuation 
spectrum depends in some way on the initial state. The 
question is how strongly. We need to understand quan¬ 
titatively how much of the Initial Information survives at 
the time when the relevant primordial fluctuations go out 
of the horizon. This general type of problem is addressed 
within non-equilibrium quantum field theory. We under¬ 
take here first steps within the conceptual framework of the 
quantum effective action. This allows us to discuss inter¬ 
acting fields and to go beyond investigations of free 
fields in a pure quantum state. Interactions are crucial for 
possible symmetrization and equilibration processes. 

The quantum effective action T for gravity, coupled to 
scalars and other fields. Is the basis for the cosmological 
field equations. Indeed, the field equations derived by vari¬ 


ation of r are exact, in contrast to those derived from the 
“classical action”, which need to be corrected by fluctua¬ 
tion effects. The issue of quantum gravity can be viewed 
as the problem of computing the quantum effective action. 
String theory or the functional renormalization group Q, 
[13 among others, have made important progress in 
this direction. 

The effective action provides essential information be¬ 
yond the exact field equations for the mean values of fields. 
It describes all aspects of fluctuations around the mean 
field as well. The second functional derivative T*^^^ of the 
effective action constitutes the exact inverse propagator. 
If it is invertible within a suitable space of field functions, 
the Inversion of T^^^ provides for a unique connected two- 
point correlation function (propagator, Green‘s function). 
This relation is exact. In a cosmological setting the equal 
time Green’s function comprises the full Information about 
the power spectrum of the cosmic fluctuations. Further¬ 
more, the third functional derivative T^^^ of the effective 
action is the one-particle irreducible three point function, 
from which the bispectrum can be extracted. Indeed, the 
connected three point function (bispectrum) obtains from 
by multiplication with three propagators and amounts 
therefore to a tree calculation. All loop effects are already 
incorporated in the computation of T. Higher order corre¬ 
lations can be obtained from higher functional derivatives 

of r. 

Once r Is computed or assumed, the propagator is deter¬ 
mined uniquely only if T*^^^ is invertible. For massless fields 
invertlbillty typically requires a regularization. This may 
be provided by a small mass term or, more generally, by an 
effective Infrared cutoff as, for example, within the effective 
average action Q. Special care Is needed In case of local 
continuous symmetries. Typically, local symmetries need 
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some gauge fixing procedure. For non-compact geometries 
or spaces with boundary invertibility needs a specification 
of boundary conditions. Once the issue of invertibility of 
is taken care of, its relation to the propagator G is an 
exact identity, 

r(2^G = i. (1) 

In this work we will assume a given form of the effective 
action and compute the propagator by solving eq. (HD- 
This equation constitutes a differential equation for the 
time evolution of the Green’s function which is solved as 
an initial value problem. All our computations are uniquely 
based on the solution of eq. o, no other input is employed. 

We emphasize that no further corrections of any kind 
are needed once F is given. All what is usually considered 
as higher order corrections in a perturbative framework, 
or quantum corrections in a semi-classical approach, is al¬ 
ready contained in the effective action. The ability of find¬ 
ing the Green‘s function with a certain precision is directly 
connected to the ability of computing F to the required 
precision. Powerful functional methods are available for 
the computation of the effective action. Often general as¬ 
pects, such as such symmetries and a derivative expansion, 
are sufficient to determine its main features. 

For a given effective action cosmology reduces to a set¬ 
ting of classical field theory. This concerns both the field 
equations for mean values and the fluctuations. We will not 
encounter any operators, commutation relations or quan¬ 
tum vacua in our discussion. On the level of the effective 
action there is no distinction between classical and quan¬ 
tum fluctuations. Only one object describes all aspects of 
Gaussian fluctuations - it is the Green’s function as ob¬ 
tained from the second functional derivative of F. This 
raises interesting questions on the “quantum origin” of the 
cosmic fluctuations. In quantum field theory F is regarded 
as a quantum object - hence the name “quantum effective 
action”. It can also be viewed, however, as a generating 
functional obtained by a Legendre transform of the loga¬ 
rithm of the partition function in a rather general setting 
of statistical physics. This poses the question if and which 
particular quantum properties are needed for the under¬ 
standing of the origin of cosmic fluctuations. 

Field relativity [il| states that the expectation values of 
physical observables do not depend on the choice of fields 
used for their description. Observations are the same for 
all “frames” obtained from each other by arbitrary non¬ 
linear field transformations. This strong principle is only 
practically valid on the level of the effective action, where 
observables are simply expressed as functionals of fields or 
correlation functions. (On the level of functional integrals 
a field transformation involves a Jacobian that is most of¬ 
ten too complicated to be handled in practice.) In order to 
establish the equivalence of different frames for fluctuation 
properties it is therefore crucial that those can be extracted 
from the effective action. This is an additional motivation 
for the present note. On the level of the effective action 
the expression of a given fluctuation observable in terms 
of fields transforms according to the general rules of field 
transformations. For example, the correlation function is 


a field bilinear and transforms as such. Typically, the cor¬ 
relation functions themselves are not frame-invariant, but 
the physical observables describing the properties of fluc¬ 
tuations are. Within the framework of the effective action 
the mapping between different frames becomes much easier 
than the question which quantum state in one frame corre¬ 
sponds to the vacuum (or other quantum state) in another 
frame. 

The standard treatment of primordial cosmic fluctua¬ 
tions [l^ - fisjl employs the linearized field equations for 
small deviations from a background cosmology. They per¬ 
mit to compute the time evolution of the power spectrum, 
which is associated to the squared fluctuations in momen¬ 
tum space. As for any linear homogeneous problem an 
initial value is needed in order to fix the amplitude. This 
is usually provided by invoking the quantum fluctuations 
in a given “vacuum” for a given cosmological geometry, as 
for example the “Bunch-Davies” vacuum [l^ for de Sitter 
space. There is an ongoing discussion on the selection of the 
“correct vacuum” for de Sitter space, [20l - [^ . More gener¬ 
ally, the question has to be answered why the fluctuations 
are described by the vacuum and not by some arbitrary ex¬ 
cited state, mixed quantum state or classical state. From a 
pure conceptual viewpoint an arbitrary consistent “initial” 
state is allowed. Properties of the fluctuation spectrum 
may depend strongly on the assumed initial state. 

Within the operator formalism the amplitude of the fluc¬ 
tuations is (partially) fixed by commutation relations for 
creation and annihilation operators. In our treatment the 
equivalent properties follow from the inversion of F^^^. The 
defining equation (HD is not homogeneous due to the non¬ 
vanishing right hand side, such that the normalization of 
G is not arbitrary. 

Within the concept of the functional integral underlying 
the effective action, the “initial values” are set as bound¬ 
ary values of fields and correlation functions. They may 
be set at the “beginning of the universe”, which actually 
corresponds to infinite physical time in the past for typical 
inflationary cosmologies [^. We are less ambitious here 
and set boundary conditions or “initial values” at the be¬ 
ginning of the inflationary epoch. For certain models this 
may coincide with the “beginning of the universe” in the 
infinite past [ll|, [1^. The issue of the influence of initial 
values can now be addressed in a practical way as the ques¬ 
tion to what extent information about the boundary values 
of correlations is still available for some given cosmological 
time. 

As a general property, such information has a tendency 
to be “forgotten” as time increases. One often finds an ap¬ 
proach to a symmetric state. The lower functional deriva¬ 
tives of F are then largely fixed by a few couplings con¬ 
sistent with symmetries and effective locality (derivative 
expansion). Assuming that for a given epoch this sim¬ 
ple form of F is already (almost) realized, the information 
on boundary conditions is no longer encoded in F*-^) on a 
practical level of accuracy. As a consequence, F*^^) is not 
invertible and the correlation function is no longer fixed 
uniquely by F^^^. In this approximation the effective ac¬ 
tion determines only the time evolution equations for the 
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correlation functions. The correlation functions themselves 
are uniquely specified only once the initial conditions for 
the time evolution are fixed. The information about the 
initial state appears now in the form of initial values for 
differential equations. This is similar to the evolution of 
the mean fields. A given T determines the field equations. 
Initial values are needed to select a cosmological solution. 
In this paper we assume an effective action that no longer 
shows a dependence on initial conditions. In particular, 
it does not explicitly depend on time. We therefore treat 
the solution of eq. © as an initial value problem for a 
differential evolution equation. 

A central part of our investigation are therefore time 
evolution equations for correlation functions. They permit 
us to start with arbitrary initial conditions for the fluctua¬ 
tions and to follow them until horizon crossing and beyond. 
In principle, the effective action and the initial conditions 
determine the correlation function for all times. Suitable 
composite fields for the energy momentum tensor of mat¬ 
ter and radiation may be needed for an efficient description 
of late cosmology. We are concerned here mainly with the 
early stage of inflationary cosmology where radiation and 
matter play no role. The effective action is then typically 
a functional of the metric and a scalar inflaton field. In the 
present paper we simplify even further by considering the 
correlation function of a scalar field in a given geometric 
background. For a given effective action for the scalar field, 
a given geometric background, and given initial values for 
the Green’s function the propagator is then fixed for all 
times. 

While the general structure of the evolution equation is 
exact, our knowledge of the effective action is at best a good 
approximation. Approximative forms can be computed us¬ 
ing appropriate non-equilibration quantum field theoreti¬ 
cal methods in the in-in-formalism [2^ Our task is 
more modest here - we will compute the consequences of a 
given assumed form of the effective action. For the effec¬ 
tive action of the scalar field we mainly consider for practi¬ 
cal computations the standard form of a covariant kinetic 
term plus a potential. These are the lowest order terms 
in a derivative expansion. In this setting we derive an ex¬ 
act time evolution equation for the Green’s function. For 
a subset of initial conditions this equation is equivalent to 
the standard approach of computing the linear time evo¬ 
lution of a small fluctuation around the background, but 
shown to be exact in our context. 

We argue that such an exact evolution equation holds 
in much more general circumstances far beyond the lowest 
order derivative expansion. Indeed, we expect the effective 
action to contain higher order terms, for example reflecting 
the scale anomaly (see ref. [l^). Such terms are typically 
subleading for our considerations since they involve higher 
order curvature invariants and higher scalar derivatives. 
Once additional local or non-local terms for the quantum 
effective action can be computed reliably, their effect can 
be included in our treatment in a well defined way. 

Within the present approximation the physical outcome 
of this paper is in many aspects equivalent to previous com¬ 
putations. There is an important practical point beyond 


previous results, however. We discuss the most general 
solution for the correlation function. It includes mixed 
quantum states in addition to the previously considered 
pure quantum states. As a result, initial states differing 
from the Bunch-Davies vacuum need not to show oscil¬ 
latory behavior. Particle or entropy production are not 
necessary. Gonceptually, the inclusion of interactions for 
the initial conditions and the general evolution of corre¬ 
lation functions allow us to discuss for a first time how 
possible equilibration and symmetrization processes could 
occur. Treatments of free quantum fields omit such effects 
from the beginning. 

The solution of the evolution equation for the simple 
case of a free massless scalar field in de Sitter space is 
shown in Fig. 1, where we plot the equal time correlation 
function k^G{k, t])/Hq as a function oiu = krj = —k/{aHQ) 
for different initial conditions. Here G is the correlation 
function, ry conformal time, k the wave number, a the scale 
factor and Hq the Hubble parameter. In our units the 
solutions do not depend on k. The solutions for different k 
are independent. 

The smooth curve in Fig. [T] is the scaling correlation 
corresponding to the Bunch-Davies vacuum [l9|. The os¬ 
cillating curve that vanishes for u —> 0 and touches zero 
at regular intervals corresponds to the Green’s function 
which is well defined in position space and consistent with 
the symmetries of de Sitter space, cf. ref. [2^. The 
two other curves employ more arbitrary initial conditions. 
The generic behavior shows oscillations whose amplitude 
decreases due to Hubble damping. At horizon crossing 
{u = —1) the relation between fluctuation amplitude and 
Hubble parameter depends on the initial conditions. 

It is obvious from Fig. I that information about initial 
conditions is not lost at the time of horizon crossing. It is 
this moment that is relevant for observations, since in the 
presence of metric fluctuations there exists a scalar fluctua¬ 
tion quantity that remains preserved after horizon crossing 
and transports the information of the primordial spectrum 
to the density fluctuations in the late universe. Our find¬ 
ing of a preserved memory of initial conditions differs from 
a tendency towards a loss of memory discussed in refs [^, 

. This is due to the facts that we discuss the fluctuation 
spectrum rather than the energy momentum tensor, that 
the relevant time is horizon crossing and not asymptotic 
time, and that we discuss initial conditions that are less 
restrictive, not imposing a pure quantum state. The ab¬ 
sence of a loss of memory is implicit in the treatment of free 
quantum fields in ref. . We show here that memory of 
initial conditions is preserved for a rather wide setting of 
models with interactions, provided that explicit boundary 
terms in F and backreaction effects for the evolution of the 
mean fields can be neglected. 

As mentioned already, the overall amplitude of G is set 
by the initial conditions independently for every /c-mode. 
Even if we restrict our attention to non-oscillatory solu¬ 
tions, important memory of initial conditions is kept by 
the time evolution of correlation functions. This concerns 
both the amplitude and the shape (spectral index) of the 
power spectrum. Indeed, the memory of initial conditions 
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has two aspects. We may not impose de Sitter symme¬ 
try on the initial conditions because we want to study a 
possible symmetrization starting from more arbitrary ini¬ 
tial conditions. In this case the evolution shown in Fig. 1 
holds for each momentum mode separately. For each k the 
overall amplitude of G can be multiplied by an arbitrary k- 
dependent factor. One finds that no symmetrization occurs 
in our approximation. In particular, the spectral index of 
the fluctuations Us receives a contribution from the shape 
of the initial spectrum. It is no longer uniquely determined 
by the shape of the inflaton potential. 

One may also study what happens if de Sitter symmetry 
or other symmetries of the mean field solution are imposed 
on the initial conditions for the correlation functions. In 
this case the amplitudes for different /c-modes are related by 
symmetry. For Minkowski space the requirement of sym¬ 
metry is sufficient to make F^^^ invertible and therefore to 
obtain a unique Green’s function. This is not the case for 
de Sitter space, mainly due to the lack of time translation 
invariance. The symmetries of de Sitter space are not suffi¬ 
cient to fix the propagator of a scalar field uniquely, even if 
we regulate by a small mass term [2^ [2l| . In the presence 
of de Sitter symmetry the general solution to the propa¬ 
gator equation involves three free parameters which reflect 
the dependence on the initial conditions. These parame¬ 
ters can influence the amplitude of the primordial fluctua¬ 
tions substantially. If nothing fixes these parameters, the 
observable cosmic fluctuations retain information about a 
particular initial state even if de Sitter symmetry of the 
correlation function is imposed. 

Different criteria for a de Sitter invariant propagator lead 
to different correlation functions at the time of horizon 
crossing. If one requires that for de Sitter space a well de¬ 
fined propagator exists in position space and is consistent 
with the symmetries, the spectrum of cosmic fluctuations 
is modified as compared to the standard scenario. While 
the shape of the spectrum (spectral index) reflects mainly 
the symmetry of de Sitter space and remains unaltered, the 
amplitude is reduced by a factor around ten. 

The absence of a loss of memory visible in Fig. [1] has 
been derived with the lowest order derivative expansion 
for the effective action. It actually holds for a much wider 
class of effective actions. The basic reason is that eq. © 
restricts the normalization of G only partially. As long as 
F does not depend explicitly on G, and the evolution of 
the mean fields is independent of G (neglected backreac- 
tions), eq. (Hj) remains linear in G. In this approximation 
the evolution equation for the equal-time correlation func¬ 
tion will turn out to be a homogeneous linear equation, 
such that the amplitude remains arbitrary. This feature is 
exact for non-interacting free fields. In the presence of in¬ 
teractions, however, terms non-linear in G will generically 
appear in the evolution equation. Beyond our approxi¬ 
mation (which neglects such terms) they may eventually 
induce “equilibration” and lead to a loss of memory of the 
detailed initial conditions. We suggest that this loss of 
memory actually happens if inflation lasts sufficiently long 
before the horizon crossing of the observable fluctuations, 
such that equilibration has enough time to occur for the 




FIG. 1: Time evolution of correlation function for different ini¬ 
tial conditions (see text). We show the equal time Green’s func¬ 
tion G{k,T]) in units of Ho and multiplied with —fc®, as a func¬ 
tion of u = —k/{aHo). In these units the curves do not depend 
on the comoving wave number k. Except for the smooth “uni¬ 
versal de Sitter propagator” the characteristic behavior shows 
Hubble-damped oscillations. Memory of initial conditions is 
kept at horizon crossing for u = — 1. This is best seen in the 
lower part of the figure which focuses on a smaller interval in u. 


observable modes. This will then select a unique form of 
the correlation function for the relevant range of scales. 

More precisely, we conjecture that for a sufficiently long 
time between the beginning of inflation and horizon cross¬ 
ing of the observable fluctuations the time evolution se¬ 
lects a unique “asymptotic propagator” or “scaling correla¬ 
tion” . Rather arbitrary initial conditions will be attracted 
towards the scaling correlation. The precise form of the 
scaling solution may not always be simple, however. 

A condition for the selection of the scaling correlation is 
formulated by analytic continuation. We postulate that the 
analytically continued Green’s function should not increase 
exponentially in the infinite past. Our condition is moti¬ 
vated by the following observations. For Minkowski space 
it selects uniquely the Lorentz invariant vacuum propaga¬ 
tor. For a long enough duration of inflation before horizon 
crossing the observable modes are far inside the horizon 
during the early stages of inflation, k/a H. Flat space 
becomes a very good approximation for these modes. Our 
conjecture holds if there is enough time such that the cor¬ 
relation function for these modes can equilibrate to the 
Lorentz-invariant propagator. Subsequently, fc/a decreases 
and the correlation function will deviate from the Lorentz 
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invariant form. We will show, however, that our condition 
concerning the analytic continuation remains preserved by 
the time evolution. It is therefore sufficient that the scal¬ 
ing correlation is reached at some early stage of inflation. 
The propagator is then given by the scaling correlation for 
all later times. For a de Sitter geometry and a derivative 
expansion of the effective action our criterion leads to the 
correlation function of the Bunch-Davies vacuum. It can be 
formulated, however, even if no exact symmetry is present. 

While the description of the equilibration process needs 
non-linearities in G beyond our approximation, the scal¬ 
ing form of the propagator can be computed for simple 
approximations of the effective action. One only needs to 
impose the condition of analytic continuation on the space 
of functions G. This makes T^^^ invertible (for proper reg¬ 
ularization of massless fields and local symmetries) such 
that eq. ([T|) has a unique solution. Indeed, if our condition 
for the scaling correlation is valid, the latter is determined 
uniquely for a given form of the effective action. For a 
derivative expansion of the effective action for the scalar 
field with up to two derivatives, a Bunch-Davies type cor¬ 
relation is selected. Modifications occur for a more compli¬ 
cated form of the effective action. In particular, this may 
concern the long wavelength tail of the correlation function. 
For the short wavelength behavior, which is relevant for the 
observable fluctuations, we expect such modifications to be 
small. 

Equilibration towards a scaling correlation is not possi¬ 
ble for free fields. The inclusion of interactions is therefore 
crucial for any possible equilibration process, (investiga¬ 
tions of the effects of interactions in a non-flat geometry 
can be found in refs. [13, In our simple approxi¬ 

mation interactions appear in the form of the inflaton po¬ 
tential. Our approach allows to study their effect not only 
for the evolution of mode functions and correlations but 
also for the initial conditions (“vacuum”). Interactions are 
per se not yet sufficient to guarantee equilibration, how¬ 
ever. As mentioned already, equilibration needs effective 
non-linearity in G for the evolution equation. 

The question if an asymptotic attractor exists, and if 
it obeys our conjecture concerning the properties of an¬ 
alytic continuation, can only be settled by an investiga¬ 
tion of the time evolution of correlation functions for a 
large space of initial conditions. In the absence of explicit 
boundary terms in F, neglecting backreaction and assum¬ 
ing only terms with up to two time derivatives we find no 
approach to such a scaling correlation that is reached be¬ 
fore or at horizon crossing. With these approximations, 
however, we would not find an approach to the Lorentz 
symmetric propagator in Minkowski space or to thermal 
equilibrium either. Even though not found yet explicitly, 
we expect some type of equilibration process that brings 
the correlation functions at least for ultra short distances 
to the Lorentz invariant propagator. For a very long dura¬ 
tion of inflation before the observable fluctuations cross the 
horizon it is the information about the initial ultra short 
wavelengths that becomes observable. A very long epoch of 
inflation can bring them to a Lorentz invariant form even 
if the equilibration time is extremely long. This will be 


sufficient to establish our conjecture at least for an infinite 
duration of inflation. 

For a given inflationary model with a finite number of 
e-foldings A^in between the beginning of inflation and hori¬ 
zon crossing of the observable modes the practical ques¬ 
tion not solved in the present paper is the determination 
of the number of e-foldings that is needed in order to 

achieve “equilibration” to the universal scaling correlation. 
In view of the small scalar self-interactions and the weak¬ 
ness of gravitational interactions we suspect that is 

rather large. For a given inflationary model is the 

crucial quantity in order to judge if the memory of ini¬ 
tial conditions is effectively lost {Nin > or if initial 

information is still observable (A^m < in the CMB- 

spectrum. 

In this paper we illustrate the considerations above by a 
computation of the scaling correlation for a scalar field in a 
cosmological background close to de Sitter space. Several 
of the concrete results are known - our emphasis lies here on 
the derivation from the effective action and on statements 
to which extent results are exact or only approximative. 
With boundary conditions specified by the scaling corre¬ 
lation, F^^) is invertible and the Green’s function follows 
from matrix inversion. We find that the equal time corre¬ 
lation function for a massless scalar obeys in momentum 
space, with scale factor a{t) and Hubble parameter iL(t), 


where 


G(fc,a) 


1 (\. if(.y) | 2 \ 

2ka'^ V y J 




( 2 ) 

( 3 ) 


The function /(y) is close to one for small v as appropriate 
for inflation. Its precise form is computed by a numerical 
solution of a nonlinear second order differential equation. 
For de Sitter space, i' — 0, one has / = 1. While for re¬ 
alistic inflation the scalar fluctuation spectrum is modified 
by the inflaton potential and mixing with gravitational de¬ 
grees of freedom, eqs. ([2]), dS]) can actually be used for the 
tensor fluctuations (up to overall normalization). 

This paper is organized as follows: In sect II we re¬ 
call briefly the relevant properties of the quantum effec¬ 
tive action. We formulate it in a framework that permits 
a straightforward analytic continuation by changing the 
value of a “background vierbein”. We turn in sect. Ill 
to the correlation function for a massless scalar field in a 
homogeneous and isotropic cosmology. The defining “prop¬ 
agator equation” relates the Green’s function G to the sec¬ 
ond functional derivative of the effective action F^^^. It is a 
partial differential equation for G depending on conformal 
time rj and comoving space-distance r. As an example, we 
solve it in Minkowski space. In appendix A we relate the 
general short distance behavior of the propagator to the 
one in flat Minkowski space and discuss the form of the 
next to leading terms in position space. 

The propagator in three-dimensional momentum space 
is addressed in sect. IV. The propagator equation can be 
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solved for each mode independently , with an inhomoge¬ 
neous term restricting partially the normalization for each 
mode separately. Within a derivative expansion of the ef¬ 
fective action we discuss the general solution of the prop¬ 
agator equation in de Sitter space. It shows no loss of 
memory. We explore the properties of a possible unique 
scaling correlation. The symmetries of de Sitter space are 
not strong enough to select a unique propagator. For se¬ 
lecting the scaling correlation we formulate our conjecture 
on the properties of the analytically continued propagator. 
It is assumed to be realized asymptotically as a result of 
the time evolution of correlation functions. It states that 
the analytically continued Green’s function is restricted not 
to diverge exponentially towards the past. On the space 
of functions with such a boundary behavior, the second 
functional derivative F^^^ becomes invertible. The de Sit¬ 
ter propagator is then unique for the modes with non-zero 
momentum and corresponds to the Bunch-Davies vacuum. 
Our conjecture is the analogue of the usual le-presentation 
in flat space to arbitrary geometries. In appendix A we dis¬ 
cuss that the Fourier transform of this propagator is not 
well defined due to the lack of regularization for the zero- 
momentum mode. The proposed unique scaling correlation 
differs from a possible de Sitter invariant propagator that 
is well defined in position space. 

Sect. |V] addresses the central formulation of a time- 
evolution equation for the correlation function that allows 
for an investigation of the role of initial conditions. We dis¬ 
play an exact evolution equation for a free massless scalar 
held in de Sitter space and discuss the properties of its 
general solution. We hnd that the memory of initial con¬ 
ditions is not lost in this case, cf. Fig. 1. We argue that 
non-linearities due to explicit boundary terms in F or back- 
reaction are needed in order to realize our conjectured scal¬ 
ing correlation. The exact evolution equation is extended 
to arbitrary homogeneous and isotropic cosmologies and 
interacting scalar helds. The approximations concern then 
only the precise form of the effective action. We discuss the 
general structure of the evolution equation for interacting 
theories in the absence of boundary and backreaction ef¬ 
fects and show that their general solution keeps memory of 
the initial conditions. In appendix B we present the gen¬ 
eral structure of the solution of the evolution equation. In 
sect. VI we formulate effective quantum helds for interact¬ 
ing theories and use them to discuss the structure of initial 
conditions. 

In sect. lVIll we compute the universal scaling correlation 
for realistic inhationary geometries, assuming that our con¬ 
jecture is true. We make contact to the standard results 
of slow roll inhation. Our results are partly exact, as for 
constant H /and H~^d^V/dip^, with V{ip) the inhaton 
potential. The infrared divergence of the de Sitter propa¬ 
gator is regularized by a mass term. In appendix D we turn 
to the correlation function for a massive scalar in de Sitter 
space. Appendix E addresses the issue of symmetries, in 
particular the consequences of the symmetry of de Sitter 
space and coordinate scale symmetry. The universal prop¬ 
agator is assumed to be consistent with the symmetries of 
the action and mean held solution. De Sitter symmetry is 


a strong constraint, but not sufficient to hx the propagator 
uniquely. 

In sect. IVIIII we turn to the impact of initial conditions 
on the spectrum of cosmic huctuations as observed in the 
cosmic microwave background. We formulate the evolution 
equation for correlation functions for an arbitrary inhaton 
potential and discuss the general solution. The initial spec¬ 
trum is only processed at horizon crossing by the scale vi¬ 
olating effects related to the inhaton potential. If inhation 
does not last too long between its beginning and the hori¬ 
zon crossing of the observable huctuations, the neglection 
of explicit boundary terms in F as well as backreaction be¬ 
comes justihed. In this case observation can look back to 
the beginning of inhation. We present an example for an 
initial spectrum that modihes amplitude and spectral in¬ 
dex for the observed huctuations. Conclusions are drawn 
in sect. m 


II. Effective action and analytic 
continuation 

The key element for the description of huctuations is 
the propagator G{x,y) = G{r],x]r]',y), where we take for 
rj, rj' conformal time. For a homogeneous background it 
can only depend on r = x — y, G(r], x; y', y) = G(r, 77 , y'), 
with Fourier transform 



Isotropy implies G{k,y,y') = G{k,y,y'),k = |fc|. The 
power spectrum of cosmic huctuations is related to the 
equal time propagator, y = y', 

G{k,y) = G{k,y,y). (5) 

The central simple observation of this note recalls that 
the second functional derivative F*-^^ of the quantum effec¬ 
tive action yields the exact inverse propagator. If F^^^ is 
invertible the propagator and therefore the power spectrum 
can be obtained by matrix inversion. If not, the condition 

F^^^G = 1 (6) 

yields a time evolution equation for the correlation function 
whose solution depends on initial conditions. Invertibility 
depends on the space of functions that are admitted for G. 

1. Effective action 

The quantum effective action arises from a functional 
integral formulation of quantum held theory, that we recall 
briehy here for an arbitrary geometrical background and 
arbitrary signature. The starting point is the partition 
function, which is dehned by a functional integral in the 
presence of sources 

Z[J] = JV^exp(^-S + 

S = [ eL[^,e^]. (7) 

J X 
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Here the generalized vector (j){x) = (}){r],x) stands collec¬ 
tively for all fluctuating fields, J{x) are the associated 
sources and we use euclidean sign conventions for the ac¬ 
tion. The geometry is described by a “background vier- 
bein” e™, with e = det(e™). This vierbein can take com¬ 
plex values. Euclidean signature is realized for real e™, 
while a Minkowski signature follows if e™ assumes purely 
imaginary values (with real e™). In particular, for flat 
euclidean space one has e™ = S^,e = 1, whereas flat 
Minkowski space obtains for e™ = (5™, e™ = e = i. 

Analytic continuation can be achieved in a simple way 
by changing the euclidean value of e™ to with fixed 

coordinates x^ = {r],x) and fixed fields 4> [i^. (For earlier 
somewhat different approaches to analytic continuation see 
refs. 0 , 13 .) Arbitrary background metrics 

= e^e^dran ( 8 ) 

can be described in this setting. Our fixed coordinates 
(? 7 ,x) span (or a suitable subspace), with time coordi¬ 
nate 77 corresponding to conformal time. The geometry of 
a homogeneous and isotropic universe is given by 

= a(77)Jr , (9) 

with 0 ( 77 ) the scale factor and 

d In CL 

adrj = dt , % = —-— = Ha, (10) 

orj 

where t and H are time and Hubble parameter of the 
Robertson-Walker metric. Metric fluctuations can be 
included into 4>, defining a fluctuating metric 

dfiu = + hfii,. ( 11 ) 

We define the functional 

W[J]=\nZ[J], (12) 

such that 

dW 

^ = ( 13 ) 

As usual, the quantum effective action r[(()] obtains by a 
Legendre transform 

= = (14) 

with J[(()] obtained by inverting eq. (I13|) . The func¬ 
tional W[J] is the generating functional for the connected 
Green’s functions for tj). The second functional derivative 
= d^r/d(j){x)d(j){y) and the second functional deriva¬ 
tive = d'^W/dJ{x)dJ(y) can be viewed as matrices. 
The matrix identity 


r( 2 )p^( 2 ) ^ ^ 


follows directly from the Legendre transformation (HI- 


2 . Correlation function 

The correlation function (connected two point function, 
Green’s function, propagator) is defined as 


G{x,y) = {^{x)^{y)) - {^(x)){^{y)) 


d^W 

dJ{x)dJ{y)' 


(16) 


(In the operator formalism this corresponds to a time or¬ 
dered operator product.) By definition the Green’s func¬ 
tion is a symmetric 


G{y,x) = G{x,y). (17) 

If r(2) is invertible the Green’s function can be found by 
matrix inversion of the “inverse propagator” 

(a<^((r)a<(.(77)) • 

In this case the propagator can be extracted directly from 
the quantum effective action, without any further assump¬ 
tions. 

The effective action is well defined for fields and sources 
for which the relation (j)[J] is invertible. The relation (1151) 
requires T^^^ to be a regular matrix without zero eigenval¬ 
ues. This is typically only the case if explicit boundary 
terms are included in T or if boundary conditions are im¬ 
posed on G. Formally, this is achieved by adding to S 
in eq. © boundary terms that are non-vanishing only at 
the boundaries. We will assume such boundary terms in 
the infinite past 77 —>■ —00 or for large negative 77 . They 
are equivalent to a specification of the initial values of 
fields and correlation functions in the infinite past [i^ . 

The boundary conditions affect the form of T close to the 
boundary. With restrictive enough boundary conditions 
f[J] is unique and T*^^) is invertible. 

In this paper we do not deal explicitly with the bound¬ 
ary terms in S or T. We rather consider a finite time 77 
for which we assume that T has already a rather simple 
form. We will question the validity of this approximation 
in sect. V. Eq. (jB]) is valid for all 77 . Without the ex¬ 
plicit boundary terms T^^^ is no longer invertible, however. 
We have therefore to consider the general solution of eq. 
([ 6 |l which involves initial values as free parameters. The 
effect of a given set of boundary terms translates to the 
selection of a given set of initial values. At this stage we 
may still impose generic boundary conditions as the ab¬ 
sence of a too strong divergence of correlation functions 
towards the infinite past. Furthermore, one may require 
that the symmetry of the background solution is shared by 
the correlation functions. 

For euclidean signature T*^^^ often becomes invertible if 
one excludes functions that diverge too strongly at the 
boundaries. This selects a unique “euclidean correlation 
function”. For Minkowski signature the justification of 
boundary conditions is less obvious. Our conjecture for 
a unique scaling correlation concerns a specification of the 
boundary conditions for the analytically continued Green’s 
function in the vicinity of Minkowski signature. The scal¬ 
ing correlation is then the analytic continuation of the 







unique euclidean correlation function. Ultimately, our con¬ 
jecture has to be justified by a study of the time evolution 
of correlation functions with arbitrary initial conditions at 
the boundary. It holds, at best, asymptotically if sufficient 
time for “equilibration” has passed since the time when ini¬ 
tial conditions are set. The choice of boundary conditions 
is an analogue to the choice of a “vacuum” in the opera¬ 
tor formalism of quantum mechanics. We will consider here 
more general initial states, including excited pure quantum 
states, mixed states and “classical states”. 

Even the specification of generic boundary conditions 
does not render invertible for all situations. Massless 
fields often induce zero eigenvalues of Obtaining a 

unique correlation function for massless fields may necessi¬ 
tate the addition of a regulator R to , such that -|-i? 
is invertible even if is not. One may then consider the 
limit i? —?> 0. The presence of a suitable infrared re^lator 
leads to the concept of the effective average action Q . 

III. Correlation function for free mass¬ 
less scalar field in homogeneous and 
isotropic cosmology 

Before turning later to more realistic inflationary models 
with an inflaton potential and interactions we investigate 
in this section the correlation function of a free massless 
scalar field in a homogeneous and isotropic geometry. We 
start with flat Minkowski space where Poincare symmetry 
can be employed to define the unique scaling correlation. 
This scaling correlation corresponds to the propagator in 
the vacuum. It transforms under scale transformations ac¬ 
cording to its dimension and contains no parameter with 
dimension of mass - hence its name. For more complicated 
geometries scale symmetry remains no longer realized. We 
conjecture that a possible scaling correlation towards which 
correlation functions are attracted for arbitrary initial con¬ 
ditions has a high momentum limit given by the universal 
free propagator. 

Time translation symmetry or Lorentz symmetry will 
no longer act as a selection principle to find the scaling 
propagator for more complicated geometries. We can, in¬ 
stead, use the properties of analytic continuation. Indeed, 
for the analytic continuation to a flat space neighboring 
Minkowski space one finds that T^^^ is invertible provided 
one imposes the generic boundary condition that the cor¬ 
relation function should not diverge exponentially towards 
the infinite past. This criterion for the scaling correlation 
needs no particular symmetry and can be generalized to a 
wide class of geometries. For this reason we construct in 
this section the universal propagator for Minkowski space 
by analytic continuation. 

In appendix A we discuss the general short distance 
behavior of the correlation function in position space. 
For arbitrary geometries it resembles closely the one for 
Minkowski space. This does not imply, however, that the 
high momentum behavior of the propagator is the one for 
Minkowski space. 


1. Propagator equation 

The identity ([51) is the central equation for the compu¬ 
tation of the propagator of a real scalar field (p{x) in a 
homogeneous cosmological background ([9]). For the partic¬ 
ular case of a free massless field one has L = d^ipdf^ip and 
the effective action reads 

r= f ea^{e~'^dnipdr,^p + dkipdk^), (19) 

J X 

with e = i for Minkowski signature and e = 1 for the 
euclidean version, e = ea"^, = a~^6ki- 

The second functional derivative T^^^ can be written as a 
differential operator 

{x, y) = (ry, f; ry', y) (20) 

= —e5{r] — r]')5^{x — y){e~^drj'a^drj' + a^Ay}. 

The “propagator equation” (I5|) takes therefore the form 

{dy + 2ndr, + e^A^)G{x,y) = — ^6{T]-r)')S^{x-y). ( 21 ) 

Our task is a solution of this equation for the Green’s func¬ 
tion G. 

The space of functions G is restricted by exchange- 
symmetry HID and “reality”. For euclidean signature a 
real action implies that G is real. For Minkowski signature 
we write 

G{r],x;r]\y) = Gsiv,x;r]',y) (22) 

+ Ga{r],x;T]',y)[9{'n-r]')-e{r]'-y)], 

with step function 9{x) = 1 for x > 0, 9{x) = 0 for x < 0 
and 0(0) = 1/2. Symmetry implies that Gs and Ga are 
symmetric and antisymmetric functions, respectively, 

Gs{y',y;v,x) = Gs{y,x;y',y), 

Ga{y',y-,v,x) = -Ga{y,x;y',y). (23) 

For Minkowski signature the action and the argument of 
the exponential function in eq. o are purely imaginary. 
The reality condition for Minkowski signature states that 
Gs is real and Ga purely imaginary, 

Gt = Gs, G: = -Ga. (24) 

For Minkowski signature solutions of the propagator 
equation m can be found if both Gs and Ga obey the 
homogeneous equation 

DGs=0 , DGa = 0, (25) 

with 

D = d^ + 2ndy - A^. (26) 

With 

D(Ga[9iy-y')-9{y'-y)]) = 2a,G,|,.,A(7y - ry') (27) 
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the solution of the inhomogeneous equation (j21[) requires with 


driGa\ri=ri' — 2a'^^ 


(28) 


f = f d^k/{2'iT f' j f — f du}j2Tr. 

t/ /c t/ t/ CO t/ 


(32) 


For Minkowski signature the “defining equation” for the 
Green’s function G reads explicitly. 

ia^DG{ri,x;ri',y) = S{r] - r]')S^{x-y), 

D = d^ + 2ndr,- A,, (29) 

such that for purely imaginary Ga both sides of eq. HMD 
are real, or both sides of eq. (l28l) purely imaginary. Eqs. 
(1^ , (1^ , (IMl) , (05|) , (177)) , (1^ generalize to a much wider 
setting than a free massless scalar field. 

The propagator is uniquely defined only if the operator 
D is invertible. If not, there exist non-trivial solutions of 
the differential equation DH — 0. Any combination G + H 
solves then eq. dUl) as well as G. For Minkowski signa¬ 
ture the propagator equation ( 1 ^^ admits a large family 
of solutions. This is obvious from eqs. 025]), OMl)- The 
inhomogeneous term only affects Ga, while the differential 
equation for Ga is a linear equation without further restric¬ 
tions. Some statement on the initial conditions for Gg is 
needed in order to specify the propagator. 

We are interested in solutions that are connected to the 
euclidean setting by analytic continuation. We therefore 
choose e™ = ^-'=)a = *(1 _ ie'ja, e > 0 , and take e —^ 0 

at the end. Eq. (12^ is thus modified 

ia^(l — ie)DG = d{ri — ri')S^{x — y), 

D = {l + 2ie){d^ + 2ndr,)-A. (30) 

Our conjecture for the scaling correlation states that for 
e > 0 the Green’s function should not increase exponen¬ 
tially for 77 —>■ — 00 . With this conjecture the only remain¬ 
ing solution of DH = Q \s H = const. This degeneracy 
with respect to a constant shift in G is characteristic for 
a massless field. To make F^^^ fully invertible a regulator 
is needed. For example, the addition of a small mass term 
makes F^^^ regular, as we will discuss later. Then the prop¬ 
agator G and therefore the power spectrum for fluctuations 
follow from the unique solution of the differential equation 
dSOl), without any further assumption about a particular 
“vacuum state”. Of course, as mentioned before, our con¬ 
jecture needs ultimately to be justified by an investigation 
of general initial conditions for the solutions of eq. 023). 

2. Propagator in flat space 

Our task is the computation of the Green’s function for 
the operator D. Without loss of generality we take rj > rj'. 
The symmetry of G under the exchange rj ^ rj', x ^ y 
can then be used to infer the result for rj < rj'. We first 
consider flat Minkowski space with constant scale factor a 
and H = 0. We write the propagator in terms of its Fourier 
transform G. 


G{r],x;ri',y) = 



For the universal scaling propagator we assume translation 
symmetry both in space and time, such that G(a;, k; w', k') 
is diagonal both in frequency and momentum. For e > 0 
F^^) is invertible and inversion of the operator D is then 
done easily in Fourier space, 

(33) 

with k = \k\ and d{k — k') = (27r)^5^(fc — k'). 

The unique scaling correlation in position space obtains 
from eq. dSU 


G{v,x;r]',y) = 




a^(a ;2 - P-h 2 few 2 ) ■ 


Performing the w-integration yields 


G{r],x-,f]',y) = 


1 


2 a 2 fc 


„-ik(ri-ri') ik{x-y) -f.k{r)-ri’) 


(34) 


• (35) 


The three dimensional fc-integral can be reduced to an in¬ 
tegral over k = \k\, with r = \x — y\, 


G{y,x-r]',y') = J dksin{kr)e-^’^^^-^'^e-^’^^^-^'^ 

0 


dTT^a^ [r^ — (1 — 2ie){ri — r?')^] 

For rj — rj' ^ 00 the correlation function vanishes expo¬ 
nentially for e > 0. We may therefore also select it by the 
criterion that it should not diverge exponentially in this 
limit (which would be the case for e < 0). 

The result for general e™ = 0 < </? < can be 

found with similar methods. 


G(? 7 ,f; rj',y) 


1 

Hr'^Ax'^Ax^ ’ 


(37) 


with = {ij ~ v'jG),Axij^ = g^j^^Ax’'. For small e > 0 
this coincides with eq. (1551) . We recognize the analytic 
continuation from flat euclidean space as well as the asso¬ 
ciated generalized Lorentz-symmetry. For Minkowski space 
the requirements of Lorentz symmetry or the analytic con¬ 
tinuation are equivalent. They both fix the propagator 
uniquely (up to a constant). 

In the limit e —)> 0 the propagator (|36l) becomes 


lim G 
£->■0 


1 j. — (?7 — 77 ')^ 

dTT^a^ £->-0 ^^2 _ jj! _ Tj'YY + 4e^(?7 


rj'Y 

(38) 


The first term corresponds to Gs in eq. (1551) . For ^ 
{rj — rj'Y' the limit e —>■ 0 can be taken without further 


(31) 
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complication. In particular, the equal time Green’s func¬ 
tion takes the simple form 

Gs{f,v)=Gs{v,x-,T],y) = (39) 

It is time-independent and reads in momentum space 


G{k) = 


I 

2a2fc' 


(40) 


For r ^ 0 the imaginary part of the propagator (1381) can 
be written in the form (| 22 )) with 


Ga{v,x-r]',y) = --^;^^{S{r-T]+r]')-S{r+r]-r]')}. (41) 

It is straightforward to verify eqs. dUl) and (ESI). 

The universal propagator in Minkowski space can be ex¬ 
tended to a large class of homogeneous and isotropic cos¬ 
mologies. Assume that we move continuously away from 
flat space by switching on a non-trivial Ty-dependence of the 
scale factor 0 ( 77 ). The propagator will move continuously 
away from the free propagator. As long as the operator D 
in eq. EHl remains invertible for e > 0 the Green’s func¬ 
tion remains unique, determined by eq. ESI). A unique 
Green’s function would imply the loss of memory of the 
initial state of the universe. Our conjectured scaling cor¬ 
relation is precisely the one that is continuously connected 
to the Lorentz invariant correlation function in flat space. 

The issue of the influence of initial conditions concerns 
the impact of boundary terms for e = 0. In our discussion 
of the propagator in Minkowski space we have not been 
very explicit on the choice of boundary conditions. The 
propagator (1381) is the only one compatible with Poincare- 
symmetry and therefore singled out. Such a symmetry 
argument is no longer available if 0(77) depends on 77 . In 
particular, time translation symmetry is lost. For e > 0 
we will require that the scaling correlation in momentum 
space does not diverge exponentially as the past boundary 
(typically for 77 —?> —00) is approached. This boundary 
condition is independent of symmetries and sufficient to 
single out a unique propagator for e > 0 . 

For e = 0, however, there will appear other solutions 
that do not diverge exponentially at the boundary. They 
can transport information about initial conditions to any 
finite time 77 . It is a dynamical question if the unique ana¬ 
lytically continued Green’s function is approached rapidly 
enough such that memory of initial conditions is effectively 
lost. We will see that for a free scalar field in a given back¬ 
ground geometry this is actually not the case. At best our 
conjecture for a unique propagator can hold in the presence 
of interactions. 


IV. Propagator in momentum space 

In momentum space one can decompose the propagator 
equation into separate ordinary linear differential equations 
for each momentum-mode. We employ this in order to con¬ 
struct the most general solution for the correlation func¬ 
tion of a free massless scalar field in de Sitter space. It 


preserves the memory of initial conditions. We also con¬ 
struct the universal scaling correlation according to our 
conjecture. It equals the one found from the Bunch-Davies 
vacuum. In appendix A we discuss the Fourier transform to 
position space. We also identify the solution which corre¬ 
sponds to a well defined propagator in position space with 
de Sitter symmetry [ 2 ^ . The latter differs from the scaling 
correlation. 

1. Green’s function in Fourier space 

Due to the translation symmetry in the af-coordinate it is 
often convenient to make a Fourier transform in the three- 
dimensional subspace. On the other hand we keep the time 
77 , for which a Fourier transform is for most cases not useful 
because of the lack of time-translation symmetry of the 
background geometry. With 

(^(77,x)= [ e*'=’"^(77,fc), (42) 

Jk 

and the condition 

<^*( 77 , fc) = 0 ( 77 ,-fc) (43) 

for real fields 4>{x), the Green’s function for a real scalar 
field in Fourier space reads 


G{r],k;r]',k') = {if{T],k)>f*{'n',k')) - {(p{r],k)){(p*{y ,k')). 

(44) 

It is the Fourier transform of G{r], x] ij', y), 

G(77,f;77',y)= / [ e^<~^^-^'y^G{y,k-y',k'). (45) 

Jk Jk' 

We use translation and rotation symmetry, 

G{r], k, 77 ', k') = G{k, 77 , T]')6{k - k'), (46) 

where G{k,T],r]') obeys eq. (|4]). 

The matrix equation (EJ can be formulated in an arbi¬ 
trary basis. In (three-dimensional) Fourier space one has, 
with e = 

T^^'>{y',k';y,k)= (47) 

< 5 ( 77 ' — ri)6{k' — k)a^{ek^ — + 2779^)}. 

The defining equation EU for the propagator reads 

a^[ek^ - |(9^ + 2779,,)]G(fc, 77 , 77 ') =5{r]-rj'). (48) 

In Fourier space the propagator equation can be solved sep¬ 
arately for each fc-mode. In particular, the inhomogeneous 
term on the r.h.s. of eq. (|48p restricts the normalization 
for each fc-mode separately.The particular geometry of de 
Sitter space is given by 

77 = --, a = ^ , e = i. (49) 

7 Ho 
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For flat space the Fourier transform of the Minkowski 
propagator Gq reads 

1 


Go(fc,r?,77') = { exp { - ek{jq - ? 7 ')}g(T? - ??') 

(50) 


2a{r])a{ri')k 
+ exp { — ek{ri' — 'r])}9{r]' — ry) 


The relation 


1 


e(fc2 - —d^) [a{r])a{r]')GQ{k, rj, 77')] = 6{r] - 77 ') (51) 


Simple solutions with constant 7 /^, 

±-+ f 

are found if the scale factor obeys 

d^a = fdj^{a^). 

This is the case for de Sitter space with e 
-l/{Hor]),f = Ho, where 


( 59 ) 


(60) 
i,a = 


fixes the normalization of Gq uniquely. For more general 
geometries the insertion of eqs. 0 , (ED) into eq. (H5)) 
yields the Fourier transform of eq. (IA.4I) 

(a^ + 2nd^ - e^k^){G{k, 77, 77 ') - Go(fc, 77, 77 ')) 

= + d^n{T]))Go{k, 77, 77 '). (52) 

2. Solution of homogeneous equation 

For solutions of the propagator equation (H51) we consider 
first the ansatz 

G{k, 77, 77 ') = {r])w+ {r]')9{T]-r]') + w+ {r])w'^ {r]')0{r]'-rj), 

(53) 

with Wk obeying the homogeneous linear differential equa¬ 
tion 

DkW^irj) = 0 , Dk = d'^ + 2'Hdn - e^/c^. (54) 

This ensures for rj ^ rj' the homogeneous equation 


ICt = 


1 iHo ^ -ikr] ^,,-1- 


+ 




e—, wt = iw^y. (61) 


The Green’s function defined by eqs. ([53]) . (| 6 T|) is precisely 
the one of the Bunch-Davies vacuum 191. 

For geometries obeying eq. (|60p the explicit form of our 
ansatz for the propagator reads for 77 > 77' 

G>ik,r],r]') = Go(fc, 77, 77 ')-P Gi(fc, 77, 77 ')-P G2(fc, 77, 77 '), 

1 1 \ Jg-ek{ri-ri') 


(62) 


Gi{k,T],y) = 
G2{k,f],T]') = - 


0(77) 0(77') J 2efc2 ’ 

j2g-gfc(7)-r)') 


2e^k^ 

where we observe the relation 


Gi(fc,77,7?') = - p d^G2{k,r],r]'). (63) 

7(0(77) - 0(77')] 


DkG{k,r],ri') = 0, (55) 

in accordance with eq. (H5|) . The ansatz dMj) is not the 
most general solution, which we will discuss later. For 
appropriate it will describe the scaling correlation. 

We explore for Wk the form 

and realize that Go { k , 77, 77 ') in eq. (15(11) obtains from G for 
77 ^( 77 ) = 0. The functions y^{rj) obey the inhomogeneous 
differential equation 

dlvt + 2(H ± ek)dr,yk ± ‘^ek'Hy^ 

= ^={dr,n + n^). (57) 

It is straightforward to verify that our ansatz solves eq. 
(|5^ for 77 77'. We may take 77 > 77' and evaluate the 

difference between l.h.s. and r.h.s. of eq. dSH), 


Dk 


Vkiv) , vUv') 




+ ykiv)ykiv')] 


-{n^{r]) + d^nirj)) 


20(77)0(77')^ 


^-ek{rj-rj') _ 


( 58 ) 


Thus eq. dSH) is indeed obeyed, and similarly for 77 < 77'. 


Using this explicit form we can now investigate if the ansatz 
(1551) . (15^ indeed obeys the defining eq. (H51) . We will see 
that this is not the case for general /. 

3. Solution of inhomogeneous equation 

For the general investigation of the inhomogeneous term 
in the propagator equation (l 48 l) we denote the propagator 
for 77 > 77' by G> and similarly use G< for 77 < 77', 

G(fc, 77,77') = G>(fc, 77, 77')6»(77 - 77') -P G<{k,r],'q')9{T]' - 77). 

( 64 ) 

The symmetry of the Green’s function implies 

G<{k,r],T]') = G:>{k,r]',T]). ( 65 ) 

For Minkowski signature the reality condition takes the 
form 


Gl.{k,'q,r]') = G>{k,'q',ri), ( 66 ) 


while for euclidean signature G> and G< are real. 

We decompose G> into its symmetric and antisymmetric 
part with respect to an exchange of the time arguments, 

G >{ k , T ], r ]') = Gs { k , r ], r ]') + Ga { k , T ], r ]') 

Gs { k , T ], r ]') = i(G>(fc, 77 , 77 ')-pG>(fc, 77 ', 77 )) 

Ga { k , r ], r ]') = ^ (G>(fc, 77,77') - G>(fc, 77', 77)). ( 67 ) 
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The propagator for arbitrary rj and rj' becomes then 

G{k,r],r]') = Gs{k,r],r]')+Ga{k,r],r]'){20{r]-r]')-l). (68) 

Applying the operator Dk, eq. (|M)) . yields 

DkG = DkGg + {DkGa){‘20{r] — rj') — l) + 2drjGaS(ri — r)'). 

(69) 

The first two terms vanish by virtue of eq. (ISll). Only 
the antisymmetric part Ga can contribute to the inhomo¬ 
geneous term in eq. (j48]). In particular, the equal time 
propagator G{k, rf) = G{k, rj, rf) only contributes to Gs and 
therefore not to the inhomogeneous part. Besides the con¬ 
ditions 

DkGs = 0 , DkGa = 0 (70) 

the Green’s function has to obey 

a,Ga(fc,r;,77')|,=v =-^- (71) 

For the particular solution (15^ the time derivative of 
the antisymmetric part of Gq already produces the correct 
coefficient of the term ~ 6{r] — rj') in eq. (1481) . For this 
solution ((62l) eq. CD) is therefore only obeyed if 


4. General solution of propagator equation in de 
Sitter space 

For de Sitter space the Green’s function (p 1 ). (ED is 
not the only solution of the propagator equation (|4^ . For 
r] > rj' the most general solution of eq. ED reads 

G>ik,r],r]') = b{k,r]')w'^ (rj) + c{k,r]')w'l^ {r]), (75) 

with 'w^{rj) given by eq. (IHTll and b{k, rj'), c{k, rj') arbitrary 
complex functions. For rj < rj' we employ G<(/c, r;, 77 ') = 
G^{k,ri',r]). Applying eq. (1551) for ij < rj' implies 

b{k,r]') = a+{k)wl{r]') + a-{k)w]:{T^'), 

c{k, r]') = P+{k)w^ ( 77 ') -b fi-{k)w~ (77'), (76) 

with arbitrary complex constants cx±{k), j3±{k). We infer 

Ga (A:, ? 7 , 77 ') = ^ [a+ (fc) - /3- (A;)] 

X [wl (77)77;+ (77') - 77;+ (77)77;" (77')], ( 77 ) 

and 

5^Ga(fc, 77, 77')|,,=V = -^(a+(fc) - l3-{k)). (78) 


dfiiGi + G2)a\T]=r]' — 0- (72) 


One has 


^ 1 

r 1 

1 ) 


4efc2 ' 

Va(» 7 ) 

0(77')/ 



(G2)a = 

(73) 


and, for rj = rj', 


drj{Gl)a — 


J 2 L a (G ) =J^ 

2efc2a ’ 2efc2- 


(74) 


The condition (17^ is therefore obeyed only for / = 0 or / = 
'H/a. The second relation is precisely realized for de Sitter 
space with / = Hq. We conclude that the Green’s function 
of the Bunch-Davies vacuum indeed obeys eq. (H51) . De 
Sitter space is the only solution of eq. ED with constant 
'H/a. We conclude that for all other solutions of eq. (l60)) 
the ansatz (ISSl) obeys the homogeneous eq. (1551) . but fails 
to be a solution of the inhomogeneous propagator equation 

(H51) . 

More generally, there are many solutions of the homoge¬ 
neous equation (1551) . For example, any given solution of eq. 
(15^ for 77 ;+ ( 77 ) can be multiplied by an arbitrary complex 
constant a+. Replacing in eq. (153)) the factors 777 ^( 77 ) by 
0 + 777 +( 77 ) still leads to a solution of eq. (I55|) . It is only 
the inhomogeneous term in eq. (l48ll that fixes the values 
of a+. 


This fixes 


a+{k)-p_{k) = l. (79) 

The general solution of the propagator equation (l48l) there¬ 
fore involves three free functions a_(fc),/3+(fc) and /3-(fc), 
with a+(fc) fixed by a+(fc) = I -|- /3_(fc). 

We next employ the reality condition (I66|) . This imposes 
the restrictions 

a-{k) = PX{k) , a\{k) = a+{k) , /3* (fc) = P-{k). (80) 

We are left with three real free functions that we denote 
by 


a{k) = c7+(fc)-b/3_(fc), 
m = p+{k)+pi{k) 

7(fc) = -7(/3+(fc)-/?;(fc)). (81) 

The most general propagator in de Sitter space reads 

Gs(fc, 77,77') = ]^a{k){w~ (77)77;+ (77') -b 77;+ (77)77;^ ( 77 ')} 
-bi/ 3 (fc){ 77 ;+ (77)77;+ (77') -b 77;-(77)777- ( 77 ')} 

+ \ l { k ){ wl { r ]) wlij ]') - 777 fc (77)77;^(77')}, ( 82 ) 

and 

Ga{k,V,v') = ( 77 ) 77 ;+( 77 ') - 77 ;+( 77 ) 77 ;-( 77 ')}. (83) 
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The equal time correlation function can be extracted 
from Gg, 

+ (84) 

{“ - &) + §t ■ 

A requirement of positivity of G(fc, ry) imposes some re¬ 
strictions on a, P and 7 - the amplitude of the oscillating 
last two terms should not exceed the term ^ a. The func¬ 
tions a{k), P{k) and 7 (fc) can be related to initial values of 
G(fc, 77 ) and its time derivatives at large negative ry. They 
therefore contain the information about the initial condi¬ 
tions, as we will discuss in more detail in sect. ED We show 
G(fc, ry) for different values of a, /3 and 7 in fig. 1 . At the 
time of horizon crossing the memory of initial conditions is 
not lost. Since initial conditions can be set independently 
for each fc-mode the general solution (IMl) retains memory 
of both the amplitude and the shape of the initial fluctua¬ 
tion spectrum. 

One may investigate a particular subset of initial condi¬ 
tions for which the correlation function is invariant under 
the symmetries of de Sitter space. We will show in ap¬ 
pendix E that the symmetries of de Sitter space require 
the functions a, /3 and 7 to be independent of k. Never¬ 
theless, we have no unique Green’s function even at this 
stage, but rather find a three-parameter family of de Sitter 
invariant correlation functions. In the limit of modes far 
outside the horizon, |fc?y| <C 1 , the equal time correlation 
approaches 

= (85) 

For a ^ P this has the same /c-dependence and therefore 
the same spectrum of the fluctuations as for the usually as¬ 
sumed Bunch-Davies vacuum. Nevertheless, the amplitude 
of the fluctuations involves an undetermined factor a — p. 

The relevant quantity for the observable cosmic fluctu¬ 
ation is the value of G{k,r]) at horizon crossing, k/a = 
Hq, kr] = —1, cf. sect. IVIll There we find 

G(fcry = -1) = Chc^, ( 86 ) 

with 

Chc = Oi — P sin(2) -I- 7 cos(2) si a — 0.427 — 0.91,5. (87) 

For a, P, 7 of the order one this is only a moderate differ¬ 
ence from the value for the Bunch-Davies vacuum Chc = 1- 
On the other hand, large values of a, P, 7 can strongly mod¬ 
ify the fluctuation amplitude. An example with /3 = 7 = 0 
is discussed in sect. IVIIII 

If de Sitter symmetry is realized the spectrum ~ k~^ 
is independent of the detailed values of a,/3 and 7 . The 
main features of the observable spectrum follow already 


from symmetry and dimensionality, despite the fact that 
memory of initial conditions remains for the general solu¬ 
tion (I84|) . For symmetric correlation functions this memory 
affects the relation between the fluctuation amplitude and 
the Hubble parameter, which in turn is given by the in- 
flaton potential. For the connection between the observed 
fluctuation amplitude and the inflaton potential It will be 
an Important task to establish if there are good reasons for 
a choice a = l,/? = 7 = 0 , orif the fluctuation amplitude 
remains to some extent an undetermined dynamical quan¬ 
tity which retains information about the “initial state” of 
the universe. In appendix A we show that a Fourier trans¬ 
form to position space only exists for a = p. In this case 
the fluctuation amplitude is suppressed as compared to the 
scaling correlation by a factor around ten. 

5. Scaling correlation 

We next formulate our conjecture for the scaling correla¬ 
tion that may be realized as an effective attractor solution 
for a large range of Initial conditions. While within our 
present approximation no effective attraction to a univer¬ 
sal propagator Is found, attractor properties are possible 
beyond our approximation, as we will discuss In sect. El If 
such an attractor exists there are several reasons to identify 
it with the scaling correlation that will be discussed next. 

For the selection of the scaling correlation we propose 
an additional condition on the propagator which will make 
It unique (except for fc = 0). It is based on the behavior 
of the analytically continued Green’s function in the infi¬ 
nite past for ry' —>■ — cx). The condition requires that G 
should not diverge exponentially In this limit. This condi¬ 
tion is the analogue of the re-prescrlption for the Green’s 
function in Minkowski space. It is our conjecture that for 
a long enough duration of inflation before horizon cross¬ 
ing of the observable fluctuations the time evolution of the 
propagator approaches this unique “de Sitter propagator” 
for rather arbitrary initial conditions. This should hold at 
least for the short-distance modes that become observable 
later. 

In the limit of flat space our conjecture selects pre¬ 
cisely the Minkowski propagator as obtained from the ie- 
prescriptlon in sect. m Furthermore, our conjecture can 
be applied independently for every A:-mode. For a long 
enough duration of Inflation before horizon crossing of the 
observable fluctuations the observable modes correspond to 
very large k, with k/a 3> H. In this range the propagator 
equation equals the one of flat space up to tiny corrections. 
As we have seen in sect. ED our conjecture on the behavior 
of the analytically continued Green’s function precisely se¬ 
lects the Lorentz-lnvariant propagator for flat space. This 
motivates our condition for the scaling correlation. For 
practical purposes It Is actually sufficient that our conjec¬ 
ture holds for the large A:-modes at the beginning of Infla¬ 
tion. As we will see In sect. El the property characterizing 
our conjecture is preserved by the time evolution. 

Even though the direct association with the flat space 
propagator yields a very good motivation for the selection 
of the scaling correlation, it remains a conjecture that the 
flat space propagator is indeed approached as a result of the 
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time evolution. We recall that the relevant modes may have 
physical momenta k/a much larger than the Planck mass 
at the beginning of inflation and could involve unknown 
“trans-Planckian” physics Si- 
Let us now apply our conjecture to the Green’s function 
in de Sitter space. We need the mode functions for 
a value of e slightly different from e = i for Minkowski 
signature. They are given by eq. (1561) for arbitrary e = 
exp(f(/j). For (p = ^ — e one has to replace for w'^ in 
eq. (j6ip the factor by and similarly for 

^ gjyen T] the contributions ^ a-{k) 

and fi-{k) in eqs. (1751) . (1751) diverge exponentially for -q' —^ 
—oo. We require that the correlation between the value 
of the scalar field in the infinite past, > —oo,x'), 

and the one at given time q, p{q, x), is not exponentially 
divergent. This implies a_ = /3_ = 0. The function /3^{k) 
remains undetermined by these considerations, such that 
our condition still allows for q > q' 

G{k, q, q) = w-^ {q)w^ {q) + /3+(fc)u;+ ( 77 )^;+ {q). ( 88 ) 

This differs from the situation in flat space where the re¬ 
quirement that G does not diverge exponentially in the 
infinite future for q ^ 00 also excludes the combination 
{q)w'l{q'). For de Sitter space one has q < 0 and the 
correlation (1881) remains finite in this range. 

The property excluding nonzero /3+ is the reality condi¬ 
tion (15^ . Indeed, 

G*{k, 77, q') = w" {q')w+ (77) -b {q')w^ iv) ( 89 ) 


practical purposes the infrared singular behavior is often 
cut off by a physical infrared cutoff connected to a typica l 
characteristic length scale for a given observation . 

The analytic continuation of the propagator (1501) reads 
explicitly 

G>{k,r],v') = W^{v)wt{v') (91) 


1 ( 

r 1 

Ho\ 

( 1 

.Ho\ 

2k ' 

Va(f?) 

ek ) 

\aiv') 



For euclidean signature the analytic continuation of de Sit¬ 
ter space is a maximally symmetric space with negative 
curvature {R = —12iLg/e^), not the sphere. The Green’s 
function for this non-compact space is given by eq. m 
with e = 1. We observe that it is real and turns negative 
for a{q)Ho/k > 1. 

While appealing by its simplicity, our condition on the 
behavior of the analytically continued Green’s function re¬ 
quires justification by the dynamics. Ultimately, one has 
to understand that dynamically a unique Green’s function 
is singled out among a much more general class of time 
evolving correlation functions. This can at best be realized 
in the form of an asymptotic attractor solution that is ap¬ 
proached rapidly enough such that the memory about the 
initial state is effectively forgotten. We will come back to 
this question in the next section and in sect. VIII. 


V. Time evolution of correlation function 


is compatible with eqs. dH]), dMl) only for /3+ = 0. For 
the most general Green’s function consistent with reality, 
as given by eqs. dMl), dMD, one infers from a_ = = 0 

that /3 = 7 = 0, while P- = 0 implies a+ = a = 1. It 
is the combination of reality and the boundary behavior of 
the analytically continued Green’s function that selects the 
unique scaling correlation for de Sitter space. This will be 
called “de Sitter propagator”. 

By virtue of our conjecture we are left (for fc 7 ^ 0) with 
a unique propagator in momentum space, given for q > q' 

by 


G{k,q,q') 


\2a{q)a{q')k 2kP \a{q) 

jj2 'I 

+ ^[exp{-ifc(?7-?7')}. 



(90) 


For k = 0 the ambiguity of adding a constant in position 
space translates to the possibility of adding to G(fc, 77 , 77 ') a 
term ^ S{k). This ambiguity can only be resolved if we reg¬ 
ulate F(^\ e.g. by adding a small mass term. The absence 
of regularization translates to the absence of a well defined 
Fourier transformation to position space, as we discuss in 
appendix A. Indeed, the term k~^ results in an infrared 
divergence of the Fourier integral. 

This infrared divergence of the de Sitter propagator 
shows also up in loop calculations [I|, IdTl - fs^ . This is of no 
direct concern here since loop effects are supposed to be 
already incorporated in the quantum effective action. For 


The key question of this paper, namely if the observ¬ 
able cosmic fluctuations keep some memory of the details 
of initial conditions or not, can only be answered by fol¬ 
lowing the time evolution of the correlation function from 
the infinite (or very remote) past to the moment when the 
wavelength of the fluctuation goes out of the horizon. We 
want to know if the propagator at horizon crossing takes 
a universal form, only dependent on the inflaton potential 
and other properties of the effective action, but indepen¬ 
dent of the details of initial conditions. Or else if there 
remains memory of the initial conditions. For this purpose 
we have to investigate large families of initial conditions 
for correlation functions and find out if they are all at¬ 
tracted to a universal propagator, and how long it takes 
until memory of the initial conditions is essentially lost. 

A central tool will be appropriate evolution equations for 
the time dependence of Green’s functions. In this section 
we derive these equations for a scalar field in a homoge¬ 
neous and isotropic cosmological background. We consider 
interacting theories with an arbitrary inflaton potential in 
sect. IVIII An exact evolution equation is found if three 
approximations are made for the general form of the ef¬ 
fective action: (i) The effective action contains no explicit 
dependence on position and time due to boundary effects, 
(ii) Backreaction effects are neglected. In principle, the 
evolution of geometry and the inflaton field could be in¬ 
fluenced by the correlation function. We assume that this 
effect is negligible, (iii) The second functional derivative of 
the effective action contains no more than two time deriva- 
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lives. We find that with these approximations the memory 
of initial conditions is still present at horizon crossing. We 
briefly discuss possible symmetrization and equilibration 
effects that may arise from modifications of the effective 
action beyond these three approximations. 

For a given background geometry the propagator equa¬ 
tion (H51) is a local differential equation, in the sense that it 
only involves properties for a given time rj. Together with 
the boundary conditions it defines the propagator uniquely 
for all times. If the Green’s function becomes universal 
and the detailed initial conditions play no role, one ex¬ 
pects that G{r],r]') depends only on the properties of F^^^ 
in a restricted time interval that comprises both ry and rj'. 
In particular, the equal time Green’s function G{k, rj) only 
needs knowledge of the background geometry and its time 
derivatives at time around 77. This crucial feature of a 
universal propagator makes the correlation functions at a 
given time 77 independent of the global form of the quantum 
effective action. Only its properties around 77 are needed. 
This fact permits the computation of fluctuations even for 
rather complex geometries for which an overall determina¬ 
tion as for Minkowski - or de Sitter space would seem very 
hard. In contrast, if memory of initial conditions is kept, 
more global features of F^^^ would be needed. The range 
of time around 77 for which F^^^ needs to be known corre¬ 
sponds roughly to the time scale on which information on 
boundary conditions is lost. 

In the quantum mechanical operator formalism one dis¬ 
tinguishes between in-in and in-out correlations. On the 
level of the quantum effective action this issue is related 
to the boundary conditions for the Green’s function. In 
principle, boundary conditions could be imposed both in 
the far past and far future. In practice, however, a given 
set of initial conditions often severely restricts the allowed 
form of the correlation function in the far future. On the 
other hand, if one remains within the allowed range the 
future boundary conditions are not expected to play any 
role for the physics in the inflationary period. Within this 
restriction the in-in and in-out situation is physically al¬ 
most equivalent. We impose boundary conditions only in 
the remote past or the beginning of inflation. Formally, 
this corresponds to an in-in situation. In case that a uni¬ 
versal form of a scaling correlation is reached the boundary 
conditions play no role altogether. 

1. Evolution of mode functions and product 
form of propagator 

Our aim is the derivation of a time evolution equation 
for the equal-time correlation function G{k,ri). A given 
solution to this equation will typically be sufficient in or¬ 
der to reconstruct the unequal time correlation function as 
well. We may take de Sitter space as an example. A par¬ 
ticular solution of the equal time correlation function fixes 
the functions a{k), (3{k) and 7 (fc) in eq. (IMl) . In turn, the 
unequal time correlation function (15^ , (1551) is fixed. 

Let us first assume that at some time 770 the correlation 
function can be written in the product form (1531) . More 
precisely, this should hold if both 77 and 77 ' are close to 
Vo = (^ + ^0/2- We will show that this product form 


remains preserved, with mode functions {rj) solutions 
of the differential equation (1551) . Initial conditions for 
and drjW^ at 770 are given by w^{r]o) and driW^^ijo). The 
same holds for w'^{r]'). By definition of the Green’s 

function constructed in this way obeys the homogeneous 
equation (|55|) for r] ^ t]' . We therefore only need to estab¬ 
lish that the normalization (ED is preserved, such that our 
ansatz solves eq. (1551) . For a correlation function of the 
form (1531) this condition can be written in the form 


{v)wt{v) - 9rjW+{r])Wf^ ( 77 ) = (92) 

We introduce rescaled functions 

Vkiv) = aiv)Wkiv), (93) 

for which the normalization reads 

dri'i^k ivW iv) - d^Vk iv)Vk iv) = -e- (94) 

The differential equation (1551) translates to 

d'^v^ = {H^ + dn'H + e'^k'^)vf. (95) 


Let us now assume that ( 77 ) are solutions of eq. (IMI) and 
compute the evolution of the normalization 

A4(?7) = dnv'^ iv)v^ ( 77 ) - ir])v^ ( 77 ), (96) 

namely 

dnMk{r]) = d^v-^ ( 77 ) 7 ;+ ( 77 ) - 9 ^ 7 ;+ ( 77 ) 7 ;^ ( 77 ). (97) 

Here the mixed first derivatives have dropped out due to 
the relative minus sign of the two terms in eq. (1961) . Using 
eq. (1551) and noting that the r.h.s contains no derivatives 
of 7 ;^ one obtains 

dr,Afkiv)=0- (98) 

The normalization does not change with time. If it is given 
by eq. dMD for ?7 = ? 7 o it will remain so for all 77 as well. 
This establishes a simple exact evolution equation for prop¬ 
agators which have the product form (1551) . It is sufficient 
to solve the evolution equation for the factors . We will 
see in sect. IVII how propagators with product form can be 
associated with effective pure quantum states. 

We have seen that the scaling correlation has the prod¬ 
uct form dSSD. Consider now a propagator which has the 
product form at time 770 , and further has the property that 
its analytic continuation for e > 0 does not diverge for 
rj' — >■ —00. Both the product form and the analyticity con¬ 
dition remain preserved for all later times 77 . For the scaling 
correlation that is uniquely defined by these conditions it 
is therefore sufficient that it is realized at some arbitrary 
time 770 . The correlation function will then be given by the 
scaling correlation for all later times 77 > 770 . This stabil¬ 
ity property is an important ingredient for any asymptotic 
solution to which other solutions may converge effectively. 
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2. Construction of scaling correlation by solution 
of evolution equation 

The conservation of the properties of the scaling correla¬ 
tion by the time evolution provides us with a powerful tool 
for the construction of the scaling correlation at some time 
77 . it is sufficient to find the scaling correlation at some 
early time, e.g. 77 —>■ — 00 . The scaling correlation for later 
times obtains by solving the evolution of wf {rf). 

Assume one has found the general solution of eq. (1551) in 
some time neighborhood of rj. This is to be understood in 
the sense that the overall geometry may be complex, but 
that in the vicinity of 77 the behavior of geometry for much 
earlier or later times can be neglected for the general so¬ 
lution to the propagator equation. We need to select from 
the general solution the particular one that obeys the 
“analyticity constraint” for the scaling correlation. In gen¬ 
eral, it is a linear combination of two basis functions. The 
correct choice of v^{rj) as linear combination of the ba¬ 
sis functions of the general solution, including the correct 
normalization, can be established by connecting the local 
geometry smoothly to a geometry with known propagator 
in the past at 770 < 77 . The precise way how this interpo¬ 
lation is done is unimportant - the result is the same for 
all interpolations. This robustness with respect to choice 
of interpolation is a result of the uniqueness of the scaling 
correlation. 

As an example, we could have obtained the correct nor¬ 
malization of the de Sitter propagator for the different mo¬ 
mentum modes along similar lines. For a given fc-mode 
the inverse propagator becomes arbitrarily close to the free 
propagator in Minkowski space for times sufficiently far 
in the past, more precisely for a <C k/H^. For a —?> 0 
(or 77 , 77 ' —>■ — 00 ) we can therefore use the propagator 
G'o(fc, ? 7 , 77 ') which takes indeed the form (I53|) with 

u±( 77 ^-oo) = ^e±^-'=T ( 99 ) 

This fixes the selection and normalization of ( 77 ) for all 
later times. The result agrees with the explicit construction 
in sect. lYl 

The exact evolution equation for the scaling correlation 
provides for a solid basis for an approach that is often used 
in practice. One solves in Fourier space the time evolution 
of a small deviation 6(p{k, rf) from the background field. 
The equal time propagator is taken as G{k, rf) ^ \5‘^{k, 77 ) p. 
The proportionality factor is then determined by choosing 
“initial conditions” by comparison with the (approxima¬ 
tive) correlation function for a quantized scalar field. This 
procedure is usually considered as an approximation. It is 
actually exact if the field equations are the ones derived 
from the quantum effective action. The functions Wk{rf) 
are indeed proportional to 5ip(k,r]), with a fc-dependent 
proportionality factor. Both Wkijf) and 6(p{k,r]) obey the 
same linear differential evolution equation. 

3. Evolution equation for equal time scaling 
correlation 

We next turn to the general evolution equation for the 
equal time Green's function for arbitrary initial conditions. 


It is formulated as a system of differential equations for 
three equal-time Green’s functions for the field (p and its 
time derivative n = drjp. We first derive here the evolution 
equation for the product form (1551) and show later that it 
actually holds for the most general propagator, imposing 
no longer the product form. 

For fixed 77 ' the evolution in 77 obeys 

{dfi+2nd^ + k^)G{k,r],r]') =0. ( 100 ) 

The time evolution equation for the equal time correlation 
function G(fc, 77 ) = G(fc, 77 , 77 ) has to take into account that 
derivatives act on both factors of 

G{k, 77 ) = wf ( 77 ) 777 ^ ( 77 ) = G^^{k, rj). (101) 
We employ the definitions 

G^^{k,r]) = {r])w+(ij) + wf {T])d^w+{r])), 

G^^ik,rf) = dj^wf:{jq)d^wl(ri), ( 102 ) 

with 

Tr{r],k) = dn(p{r],k). (103) 

They describe connected correlation functions of the type 

i((7t(77, k)if*{r], k'))a + {(fiv, k)n*{r], k'))c) 

= G^^(fc, 77)5(fc - fc'), (104) 

( 77 ( 77 , k)Tt*{ri, k'))c = G^^{k, ri)6{k - k'). 

The evolution equation for the equal time correlation 
function takes the form 

drjGipipifk^ rfj — 2GT^ip(^kj Tj)^ 

dnG^^{k,r]) = G^^(fc, 77 ) - 2 'HG^^(fc, 77 ) - fc^G^^(fc, 77 ), 
d-qGT,-^{k,rf) = - 4 'HG^ 7 r(fc, 77 ) - 2fc^G^,^(fc, 77 ). (105) 

These equations bare certain similarities with the ones in¬ 
vestigated in ref. (see also refs. BUI) , which are 
based on exact evolution equations for the equal time effec¬ 
tive action [H^, There are two important differences, 
however. The present equations are based on the quantum 
effective action, which already includes the effects of quan¬ 
tum fluctuations. This explains why these equations can 
be exact for a simple given form of F, as in our case. In con¬ 
trast, the equations in ref. are based on the classical 
action and therefore comprise additional terms involving 
higher order vertices. A second characteristic feature is 
the Hubble damping related to 'H > 0, not present for in¬ 
vestigations in Minkowski space. We will see below that 
eq. (11051) does not need the particular product form (llOip 
but rather follows directly from the definition of the equal 
time Green’s function, together with eq. (11041) . 

4. Time dependence of equal time Green’s function 
in de Sitter space 

The system of equations (11051) permits us to address our 
question to what extent equal time correlation functions 
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keep memory of initial conditions. We can specify initial 
conditions for and Gtttt at some large negative tjq 

and follow the evolution towards a later time ry, for example 
at horizon crossing. In this section we will do so for de 
Sitter space with H = —1/?7, while we extend the discussion 
to more general geometries and the inclusion of an inflaton 
potential in sect.VIII. 

The combination 

= GipipGTTT^ — G^^ (106) 

can be used to construct a conserved quantity. (We use the 
symbol ac for reasons of consistency with earlier work on 
the time evolution of Green’s functions. It should not be 
confounded with the parameter a{k) in the general propa¬ 
gator in de Sitter space (15^ 1. The combination evolves 
according to 

= -maQ^ = (107) 

with solution 



The quantity is conserved. One may verify this ex¬ 

plicitly for Gip^ = w'^{v)w'l[r]) with given by eq. (ICTI) . 
With 

G — ^ -I- ^0 c — r< — ^ nOQl 

~ 2a^k ^ 2fc3 ’ “ 2ak ’ " 2 a 2 ’ 

one has 

( 110 ) 

We note that aG{k)a'^ is conserved for every fc-mode. The 
existence of an infinite number of conserved quantities is a 
clear sign that memory is not lost. It prevents the approach 
to a scaling correlation. 

In appendix B we analyze the behavior of the general 
solution in terms of a partial fixed point (scaling solution) 
and deviations from it. In this form the analysis can be 
generalized to a much wider context. It also highlights the 
role of the scaling correlation as a fixed point. If this fixed 
point is approached faster than by Hubble damping the 
memory of initial conditions would be lost asymptotically. 

For our present setting we actually know already the 
most general form of the solution to the evolution equation 
(11051) : it is given by eq. (|M)) . Indeed, it is straightforward 
to verify that eq. (I84|) solves the evolution equation. The 
three free integration constants a, /3 and 7 correspond to 
the three initial values for G^ip, Gt^ip and Gtttt- The asymp¬ 
totic value (IB.201) reflects eq. (1551) . We conclude that for 
a free massless scalar field in a given geometrical back¬ 
ground the information contained in the initial state is not 
lost. This information influences strongly the amplitude 
and spectrum of the fluctuations. The situation remains 
the same for a massive scalar field. This is similar to the 
behavior of the correlation function for a free scalar field 
in flat space. Infinitely many conserved quantities acik) 


prevent a loss of memory [s^ . In order to see equilibra¬ 
tion (e.g. the approach to a thermal equilibrium state) one 
needs approximations for which acik) are no longer con¬ 
served. For de Sitter space these conserved quantities are 
replaced by aQ^a^. 

5. Equilibration in de Sitter space 

In flat space, the situation changes profoundly in the 
presence of nonlinearities induced by scattering. For an 
interacting scalar field the correlation functions reach for 
asymptotic times the values of thermal equilibrium, pro¬ 
vided certain conditions for the energy distribution in the 
initial state are met . For intermediate times an approx¬ 
imate partial fixed point may be reached , correspond¬ 
ing to “prethermalization” [5a, It seems plausible, even 
though not proven at present, that in the presence of inter¬ 
actions a similar phenomenon of equilibration could take 
place in de Sitter space. This would substantiate the pic¬ 
ture of de Sitter space as a type of a “thermal state” with 
a “de Sitter temperature”. We note in this context that in¬ 
teractions are not restricted to scalar self-interactions. In 
the presence of metric fluctuations every model contains 
gravitational interactions. 

A possible equilibration process contains different facets. 
As a first aspect one would expect that equilibration leads 
to a state consistent with the symmetries of de Sitter space. 
In this case the functions a{k), j3{k) and 7 (fc) would all be¬ 
come independent of k. For example, the high momentum 
modes may equilibrate at early times to a scale invariant 
spectrum, which implies de Sitter symmetry (cf. appendix 
E). As a consequence, a (partially) equilibrated state with 
de Sitter symmetry would lead for later times, aH ^ fc, to 
Gpp ^ k~^, except for very special initial conditions that 
lead to f3 = a. The spectrum for k/aH 0 takes now 
a universal form, loosing the memory of the initial condi¬ 
tions. This also holds at horizon crossing where symmetry 
and dimensional analysis lead to a spectrum G ^ 
cf. eqs. (1551) . (1571) . The shape of the spectrum, in partic¬ 
ular the spectral index of the observable primordial fluc¬ 
tuations, is determined in this case by the background ge¬ 
ometry at the time of horizon crossing. In contrast, the 
fluctuation amplitude is not fixed by symmetry properties. 
It depends on the parameters a, /3 and 7 . An approach to 
de Sitter symmetry therefore looses only part of the initial 
information. 

As a second aspect one may expect that an equilibrated 
state obeys some “fluctuation-dissipation relation”. This 
links the amplitude of the symmetric part of the Green’s 
function, as given by a, j3 and 7 , to the one of the antisym¬ 
metric part which is fixed by the inhomogeneous equation. 
It is this type of relation that may finally fix the fluctu¬ 
ation amplitude. If equilibration leads to a — /3 = 1 the 
amplitude takes a universal form for k/{aH) <C 1, while 
not necessarily being universal at horizon crossing. If, in 
addition, /3 = 0,7 = 0 , the scaling solution corresponding 
to the Bunch-Davies vacuum will be reached. We briefly 
discuss the fluctuation-dissipation relation in appendix C. 

We conclude that equilibration in an interacting theory 
could lead to a loss of memory concerning the precise ini- 
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tial conditions. It is this phenomenon that is supposed to 
select the unique de Sitter propagator (ITOl) . as defined by 
requirements for the analytic continuation. What remains 
to be settled is an answer to the question if an equilibra¬ 
tion happens, and on what time scales it is realized. Only 
after answering these crucial questions one can asses the 
questions raised in the introduction, namely if the cosmic 
fluctuation spectrum and amplitude is independent or not 
of the detailed initial conditions. 

A key ingredient for the memory of initial conditions is 
the linearity of the system of equations (I105|) . We will see 
below that this linearity is maintained for an interacting 
scalar field, provided our assumptions on the absence of 
explicit boundary contributions to F and the neglection 
of backreaction hold. A possible equilibration process has 
therefore to go beyond these approximations. 

6. General structure of evolution eqnation 

For getting further insight into the initial value problem 
we should consider the quantum effective action in the pres¬ 
ence of boundary conditions, imposed for ry —>■ —oo or for 
some very large negative rjin. Not too far from the bound¬ 
ary the effective action F is not supposed to have the simple 
universal form m - some information about the boundary 
is expected to be still present. For a different form of F the 
flow equation for the equal time propagator will no longer 
have the simple form (11051) . We want to understand the 
structure of these modifications. 

For this purpose we will consider an approximation 
where higher than second derivatives with respect to ry can 
be neglected in F^^^ such that (omitting the unit matrix) 

f( 2 ) = AD^, 

Dn = d^ + 2Cdn + B, ( 111 ) 

with A, B and C arbitrary real quantities that are com¬ 
putable for given rj and k. (The function A should have 
no zeros such that XjA exists. Generalizations beyond the 
homogeneous and isotropic setting are possible.) In partic¬ 
ular, the functions A, B and C may depend themselves on 
other fields or on , G-^^p and G^^-tt , such that the system 
of equations becomes non-linear. They may also involve the 
antisymmetric part Go. , enforcing a normalization of Gg by 
the properties of Go- We next derive evolution equations 
for the equal time Green’s function for a general effective 
action that has a second functional derivative (HIIJ. 

We use the general form (IMl) and concentrate on the 
symmetric part Gg (IS71) . It obeys 

^^G,(7y,7y')=0, ^^.G,(7y,ry')=0, (112) 

where we omit the label k. For the second equation (I112p 
the differential operator D acts on rj'. It follows from the 
first equation by using the symmetry of the Green’s func¬ 
tion. The time evolution of Gg involves only the functions 
B and C. The function A appears in the normalization of 
the antisymmetric part Go- We may assume a normaliza¬ 
tion of the scalar field for which A = ia^{ri), cf. eq. (1^ . 
but this is not crucial. 


The general definition of G,^,^, Gt^^ and Gtttt is given by 
eqs. (HU), dH]) and (I104p , with G^p^p = Gs{r],T]). The 
relation 


drjGpp = 2G^p (113) 

follows directly from the definition. Similarly, one obtains 

dr^G-^p = Gt^-^ + + (pd'^(p*)c ( 114 ) 

= G^^ -b i(52Gs(?y,?y') +52,Gs(?y,ry'))|^,^^. 

Employing eq. (11121) and the definition (11111) this yields 
d^G^p = G^^-{C(?y)9^G,(7y,7yO+C(ry')5^,G,(7y,7y') 

+ \B{r])Gg{'q,r]') + is(77')G^(ry, 7y')}|^,^^(115) 

or 


d-qG-i^p — Gtttt ~ 2CGt^p — BGpp. 


A similar line of argument gives 


drjGTTTr — —'XCGTr-TT ~ 2SG,rvJ- 


(116) 

(117) 


Eqs. (11131) . (11161) . (11171) form for any given B and C in 
eq. dm]) a system of exact evolution equations for the 
equal time correlation functions. The ansatz (11171) covers 
a rather wide family of situations for interacting scalars. 

Eor C = % depending only on rj and B — k'^ we recover 
our previous setting (11051) . Eq. (11051) is therefore not linked 
to the particular product form of the propagator (I101|) , but 
holds for the most general propagator. 

The presence of interactions does not necessarily in¬ 
duce non-linearities in the evolution equation. We may 
include an arbitrary potential V{ip) for the inflaton. The 
second derivative dVjdip^, evaluated for the background 
field ip{r]), amounts to a mass term m^{r]) that may de¬ 
pend on r] through the background field. This results in 
B{f]) = -b m'^ {r])a‘^ {rj), but leaves the linear structure 
of the evolution equation unchanged. We will discuss the 
correlation function of massive scalar fields in sections I VIII 
and IVIIIl 

The general form of the evolution equations for the equal 
time Green’s function is given by eqs. (11131) . (11161) and 
(11171) . If higher derivative terms play no role the physics of 
a possible equilibration towards the unique de Sitter prop¬ 
agator is encoded in the two functions C{k, rj) and B{k, rj). 
There are several possible sources for non-linear effects in 
the evolution equation. Eirst, F*^^) typically depends on 
other fields. The evolution of these fields may depend on 
the value of G. An example is the energy momentum tensor 
associated to the scalar field fluctuations which involves G. 
It influences the evolution of the metric, inducing thereby 
an effective non-linearity in the system of equations. This 
effect is called “backreaction”. 

Second, the form of the effective action may be time- 
dependent itself as a result of explicit boundary contribu¬ 
tions to F. The functions B and C could depend directly 
on G{k, rj), or more generally on the unequal time Green’s 
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function G{k,r],rj'). Furthermore, higher order vertices 
may become effectively functions of time, and this results in 
a time dependent second functional derivative F^^^ through 
a non-trivial background field (p{ri,x). Away from equilib¬ 
rium F^^) could have an imaginary part. A third effect is a 
richer variety of background solutions. Even if F*^^^ remains 
unchanged, initial values may be such that the background 
solution is no longer homogeneous and isotropic. Since 
F^^) results from an expansion around the background the 
functions C and B reflect such a change of the background. 

We note that the formal way how non-equilibrium ef¬ 
fects modify the time-evolution of the propagator depends 
on the choice of fields. For example, the effect of radiation 
and matter in gravity may appear as an explicit source 
term in the background equations for the metric, in the 
form of the energy momentum tensor in the Einstein equa¬ 
tion. Alternatively, one could introduce a composite field 
for the energy momentum tensor, which would then be part 
of the solution for the background field. The situation is 
similar for the correlation functions. The non-linearities 
are most easily visible if and G,rir appear them¬ 

selves as fields. This is realized in the two-particle irre¬ 
ducible formalism [dol - l^ or for the equal time effective 
action . In the two-particle irreducible formalism 

it is the whole propagator G(a;, y) that appears as a field. 
Non-linearities are directly visible in the effective action 
^ 2 Pi['^,G]. The quantum effective action discussed in the 
present paper obtains from T 2 P 1 by inserting a solution 
G[ip\ of the field equation dT 2 Pi/dG = 0. 

For general C and B the evolution of obeys 

dr,{aQ^a^) = 4(-H - C){a-^?a^). (118) 

An approach to the unique de Sitter propagator requires 
that C differs from "H, except for a zero if = 1/4. The 

corresponding partial fixed point has to be attractive as y 
increases. A non-vanishing r.h.s. of eq. (11181) is necessary 
in order to avoid conserved quantities which would be an 
obstacle to equilibration. The quantity C — H stands for 
an effective irreversibility beyond the purely geometric one 
induced by the Hubble damping. A fixed point for this 
quantity requires non-linearity of the evolution equations 
and is part of the fluctuation-dissipation relation. 

7. General structure of solution for propagator 
equation 

Let us omit possible non-linearities and consider B and C 
as given functions of ij and fc, with A = ia?'{r]). This should 
be a reasonable approximation on time scales shorter than 
the (unknown) equilibration time associated to the non- 
linearities. We can construct the general solution of the 
propagator equation for the unequal time correlation func¬ 
tion for F^^) given by eq. (Hill) . Since Ga is fixed by the 
inhomogeneous term this amounts mainly to a solution of 
eq. H12I) . Since the evolution equation for the equal time 
correlation function follows from eq. H12I) , the general so¬ 
lution of the propagator equation automatically solves the 
evolution equations H13I) . H16I) . H17I) for the equal time 
correlation functions. 


Following the line of arguments in eqs. (113, dlS we can 
write the general solution in the form 

G>{k,T],r]') = a+{k)w'^ {r]') + l3-{k)w^ {r])w'f: {r]') 

+(3+{k)w'l {r])w^ [r]') + a-(k)w~^ {y), (119) 

with {rj) two basis functions from which the general so¬ 
lution to the homogeneous equation = 0 can be 

constructed, 

i’kiv) = C+W+{y) + c-w-^: {rj) , Dr,w^{r])=0. ( 120 ) 

Here w^{r]) need not yet to be complex conjugate to each 
other. The inhomogeneous equation constrains the combi¬ 
nation a+(fc) — /3_(fc), 

= ^(a+(fc) -/3_(fc))(w-(r7)u;+(7?') 



We can choose for each k complex multiplicative con¬ 
stants for such that these functions are normalized ac¬ 
cording to 

9v{wkiv)wtiv')-wtiv)Wk{v')\v=ri' 

With this normalization the inhomogeneous equation re¬ 
quires 

a+ik) - P-ik) = 1. (123) 

If w'^ and (wj!")* are independent we can choose w'j^ = 
(re/")*. If not, we can always choose linear combinations of 
the basis functions such that w'j^ = (re/")*. This is compat¬ 
ible with a purely imaginary l.h.s. of eq. H22I1 . With this 
choice of mode functions the reality condition implies 
the constraint (l80l) . The general solution to the propagator 
equation is therefore given by eqs. (15^ . (|55)) - 

Similar general solutions arise from the discussion of free 
quantum fields in a geometrical background. The function 
tpkiil) plays there the role of a wave function. For the op¬ 
erator formalism to be efficient it is crucial that the quan¬ 
tum fields are non-interacting. Eigenstates of a quantum 
Hamiltonian for interacting theories are very complicated 
multi-particle objects and a perturbative expansion around 
the free field theory becomes rapidly rather involved. In 
contrast, our formalism covers directly interacting theories 
as well. Interaction effects are accounted for by the com¬ 
putation of the effective action F. The latter can have a 
rather complicated form, involving arbitrary potentials or 
other non-linearities. Nevertheless, we have found a struc¬ 
ture analogous to free quantum fields, based only on the 
assumption that non-linearities from backreaction and ex¬ 
plicit boundary terms are absent and that higher order 
derivatives with respect to rj play no role for the second 
functional derivative F^^^. 

Our formalism covers a much wider range of physical sit¬ 
uations than free quantum fields in a given background. In 
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particular, it includes mixed quantum states. Moreover, 
our evolution equations for the correlation functions are 
exact once or, equivalently, the functions A,B and C, 
are known or assumed. The mode functions w^{rf) ap¬ 
pear in our treatment just as a convenient way to solve 
the propagator equation, not as fundamental objects. Our 
setting can be realized in a rather wide context of statisti¬ 
cal systems as discussed in sect. [iTl It applies to euclidean 
signature as well if we replace the factor —i in eqs. ( 11211 ) . 
( 11221 ) by —1 and adapt the reality condition. 

VI. Effective quantum fields and initial 
conditions 

The analogy to a free quantum field allows us to interpret 
i/'fe {rf) as an effective quantum field even in the presence of 
interactions. (This resembles the use of a free Fermi gas for 
electrons in a solid despite the presence of strong interac¬ 
tions.) This interpretation requires certain restrictions on 
the functions a{k),l3{k) and 7 (fc), however. If we want to 
associate G{k, ij, -q') with the propagator for a pure state of 
a free effective quantum field ipkiv) it to obey 

G>ik,q,q') = ^pkivWkW)- (124) 

With eqs. (11191) . (11201) this implies (for all k) a+ = |c_p, 
/3_ = |c+p, /3+ = c+cl, a- = c_c(j_. The reality condition 
is obeyed and eq. (11231) implies 

|c_p = l+|c+p. (125) 

For a, f3 and 7 one finds the “pure state conditions” 

a = |c_p + |c+p = l + 2 |c+p, 

P = c+c*_ + c)j_c_, 7 = —i (c+cd — c)j_c_) . (126) 

This implies, in particular, a ^ 1. Other constraints apply 
for P and 7 , e.g. 

/32 +^2 ^ _ 1. (127) 

For example, the case a = /3 = 1, which leads to a 
well defined de Sitter invariant Green’s function in po¬ 
sition space, cf. ea. (IA.21l) . is not compatible with the 
pure state condition. For pure state initial conditions all 
deviations from the scaling solutions are oscillatory, e.g. 
a > 1 implies -I- 7 ^ > 0. For de Sitter space the pure 
states correspond to the so called a-vacua [13, Uli with 
|c_p = cosh^ a, |c+p = sinh^ a. 

There is, however, no need to impose pure state initial 
conditions. Remaining within the language of quantum 
mechanics we may explore “mixed states” where 

G>{k,q,q') = (??'))* , (128) 

i 

with Pi > 0, Si Pi = 1- This yields 

a+ = ^Pi|c(!^P , P+ = Y^PiC'+ , (129) 


and similar for a- and /3_. One has again the restrictions 
a+ > 0 , /3_ > 0 which imply that a is bounded from below 
by 1 , 

a = l + 2J2p^\cf\''■ (130) 

i 

The relation for P and 7 becomes now 

+ 7^ = . (131) 

It is easy to see that for arbitrarily large a one can realize 
/3 = 7 = 0. The correlation function needs no longer to 
be oscillatory for a > 1. This shows that deviations of 
the initial spectrum from the Bunch-Davies vacuum needs 
not to induce oscillatory features in the power spectrum 
(as assumed in ref. ii) Modifications of the initial 
spectrum are not necessarily related to particle production 
or non-adiabatic processes. For mixed states one finds an 
upper bound for /3^ + 7 ^, such that the restrictions for a 
“mixed state initial condition” are 

ct{k) > 1 , /3^(fc) -l- 7 ^(fc) < a^{k) — 1. (132) 

This still excludes states with a = /3 ,7 = 0. 

The relations (|132l) imply a minimal value for the equal 
time correlation function for each fc-mode. We demonstrate 
this for de Sitter space and 7 (fc) = 0, where the local min¬ 
ima of oscillations in eq. (IMl) occur for 

Gmin{k,q) = {a{k) — P{k)) ■ (133) 

For a given a(fc) the minimal value of Gmin is realized 
for maximal P{k) = ^a^(k) — 1, corresponding to a pure 
state. The overall amplitude a — P = a — \/a^ — 1 is larger 
than one for a > 1. The lower bound on G{k,q) reflects 
the Heisenberg uncertainty relations associated to quantum 
systems - they imply a non-zero minimum of fluctuations. 
Among the pure states the one with minimal fluctuations 
is the Bunch-Davies vacuum with a = l,/3 = 7 = 0. Thus 
our conjectured scaling correlation is also singled out as an 
effective quantum vacuum with minimal fluctuations. This 
holds even in the presence of interactions. 

It is not clear if for a general interacting theory the con¬ 
straints (11321) remain necessary. Formally, their violation 
only means that no description of the correlation function 
in terms of an effective free quantum field is possible. For 
a weakly interacting quantum field theory the “fluctuation 
constraint” (11321) should hold at least approximately. We 
will in the following assume the validity of the fluctua¬ 
tion constraint. The de Sitter invariant Green’s function 
in position space (/3 = a, 7 = 0 ) can then not be realized 
exactly. For large a{k), that we may associate with high 
occupation numbers of modes, the role of the “quantum un¬ 
certainty” becomes negligible and the constraint becomes 
effectively -|- 7 ^ < . This amounts simply to the con¬ 

straint G{k, 77 ) > 0. 
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VII. Scalar correlation function in 
inflationary cosmology 


In this section we investigate the scaling correlation ac¬ 
cording to our conjecture for realistic inflationary models. 
The propagator in momentum space takes then the prod¬ 
uct from (1551) . For rj > rj' it is expressed in terms of 
(l93l) as 

G> {k, 7 ?, Tj') = {v)vt iv')- (134) 

In this form it is straightforward to compute the propagator 
for geometries different from de Sitter space. It is sufficient 
that for large negative rj and 77 ' the flat space propagator Go 
becomes a good approximation. We can then solve the evo¬ 
lution equation for ( 77 ) taking their values for flat space 
as initial condition. Furthermore, once the mode functions 
8 '!'® computed it is also straightfor¬ 
ward to construct the general solution (|55]) . (|551) . For an 
arbitrary homogeneous and isotropic cosmology it differs 
from the one for de Sitter space only by the difference in 
the mode functions. 


1. Green’s function for geometries close to de Sitter 
space 

Our starting point for the computation of the Green’s 
function for a massless scalar field in a large class of infla¬ 
tion models is the evolution equation (IMl) for v^{r]) with 
= —1. This is a real equation such that both v and v* 
are solutions. We make the ansatz 

(^) = ^^(2/)®”*^'' ’ (^) = (^fc (^))*’ (135) 


with 



(136) 


Typically, "^( 77 ) vanishes for 77—>-—00 or a—>-0 such that 
in this limit we have the propagator in Minkowski space. 
This fixes the normalization 


lim b{y) = 1 . 

y—^oo 


(137) 


The evolution equation (1551) reads for b{y) 


l + iy 


2i- 


dv 

dy 


dv - 


2 -\- V 


(1 -I- 7^)2j/2 


b{y) = 0 , 


with 




(138) 


(139) 


Here v is a slowly time dependent function defined by 

dr^n = (I + v)'H^. (140) 

For the particular case of powerlaw inflation with 
a = = aot~i 


one has a constant 1 / < 0, while de Sitter space corresponds 
to 7 ^ = 0. For more general models of inflation is a small 
quantity and obeys in a Robertson Walker metric 


v = 


H 

IP' 


In the Robertson-Walker metric one has 



(142) 


(143) 


such that modes far outside the horizon correspond to 7 / —>■ 
0 and horizon crossing occurs for y = 1. For constant v the 
asymptotic behavior for 7 / —>■ 0 is easily found by neglecting 
in eq. (I138p the term ~ namely 


\\m b{y) = bt)y i+-. 

y-70 


(144) 


This solution is an attractor for 77 —> 0, with general solu¬ 
tion 


lim b{y) = boy 1 +- ( 1 - —^y -f •.. ) -t- 61771 +- , (145) 
y-s-o 1-1-7^ / 

and 61 the second integration constant. Inserting the 
asymptotic solution (1 1441) into eqs. (IMl) . (jl35l) . (I55|) yields 
the correlation function for modes far outside the horizon 

\bo\^H^ f k 

2 \aH J 

. 

For the amplitude we first note that bg is a dimensionless 
number that only depends on ly, with bo = i for 7 ^ = 0. 
Indeed, eq. (|138ll contains no scale if du/dr] = 0. The 
constant bg obtains from the solution of this equation with 
the normalization (I137p . We may employ the ansatz 

+ (147) 

which becomes exact for 77 = 0 with f{y) = I. This yields 
the equal time propagator m mentioned in the introduc¬ 
tion. In the following we will discuss the properties of the 
function f{y) in more detail. 

The correct asymptotic behavior of b{y) both for 77 —^ 0 
and 77 —>■ 00 obtains if f{y) reaches constant values in these 
limits, with 

bo = [if{0)]^. (148) 

The differential equation for f{y) reads 

[(/ - *y)(i - iy) - iyv] dyf + ^-{dyff 

-^^^y{f - w)dyf + |(/ - iy)if - ydyf )dyiy 

= */(l-/)+2/(l-/+0- (149) 

By multiplication with i this becomes a real differential 
equation for the variable iy. For 77 —?> 0 eg. (11491) simplifies 
to 




(141) 


dyf = *(1 - /) 


(150) 
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with solution 


/ = l + ce-*^. (151) 

This tends for j/ —>• 0 indeed to a constant, 

\bo\^ = \l + c\^. (152) 

The other limit for y —>■ oo reaches for constant v 

The only solution that is consistent with the normalization 
b{y —>■ oo) = 1 occurs for 

f{y ^ oo) = 1 + ^. (154) 

Thus both 6 (y) and dyb{y) are fixed for y ^ oo and we 
can integrate eq. (I138|) unambiguously from some “initial 
value” for very large y to y —>■ 0. We find that (/ — \')lv is 
almost independent of as shown in Figl^l One obtains 

/(y = 0) = 1 + 0.273i/, (155) 

such that (for small v) 

\bo\^ = {fiy = 0))^-1 + 0.55U. (156) 



FIG. 2: Propagator near de Sitter space. We show (/ — 1)/;/ 
as function of a: = iy . The two curves for i/ = 0.01 and v = 0.5 
are almost identical. 

With eq. (11541) we can develop for y —>• oo eq. (I149|) 
iteratively, the next step being 

On the other end, for y —>■ 0, the solution (jl45l) implies 

/ = + { 1 - bl+ny = fo + - fo), (158) 

where fo = —ib]f^'^ is real for constant u. For realistic 
inflationary models \v\ is small and /(y) stays close to one. 
One therefore finds only a small correction, also in the case 
where v depends on rj, 

|6oP = l + 0(j.). (159) 

If the y-dependence of ly remains moderate eq. (11561) is a 
good guide. 


2. Scalar Green’s function far outside the horizon 

Eq. (IMI) is valid for all times (as long as our assumptions 
hold and no other degrees of freedom play a role). We 
can therefore use its solution to compute the asymptotic 
behavior of the correlation function for a massless scalar for 
infinite time - or at least times long after horizon crossing. 
The case of massive scalars will be treated in appendix D. 

For geometries close to de Sitter space the time evolu¬ 
tion for modes far outside the horizon is determined by the 
factor 


y(y) = = a-^n^ (160) 

in eq. (I146|) . Inserting for constant v eg. (11401) yields 

dy\ng{r]) = -2% + ^ = 0 . (161) 

1 -I- H 

We recover the known result that the fluctuation amplitude 
for modes far outside the horizon does not change anymore. 
To a good approximation we can evaluate the amplitude 
of the power spectrum at the time when a mode leaves the 
horizon, y = 1, where eq. (11461) yields 

k^G(k, a^oo) = = 1). (162) 


Here 


^ g{a -t oo) 
g{a = k/H) 


(163) 


is a fc-independent factor close to one. 

Corrections to this simple formula can be evaluated in 
our setting in an unambiguous way. For constant v the 
relations (11461) and (11621) are exact since eg. (11381) has no 
explicit ^-dependence and the asymptotic behavior (I144p 
fixes the y-dependence. For a varying v{rf) an additional 
fc-dependence arises from dyV = —dyiy/[{l -I- iy)k]. 

The relation (I162|) shows no explicit dependence on k 
anymore, and it depends explicitly on v only via the factor 
\bo\^G. The A:-dependence of the combination k^G{k,r]) 
arises from the A:-dependence of the value of the Hubble 
parameter at horizon crossing, H{y = 1). This is reflected 
by the shape of the pure scalar correlation. 


Tlip 1 


cHn(fc^G(fc, y)) d\nH{y=l) 
d\nk d\nk 

g dlnH dInH/dt 

d\n{aH)^y^^ d\n{aH)/dt' 


(164) 

k. 


Using eq. (11421) we find 


- 1 = T-—> (165) 

1 + iy 

in accordance with the explicit fc-dependence in eq. (11461) . 
This formula is exact for constant v. 

For realistic inflationary models eqs. (^51) . (jl35l) . (I138p 
cannot be used directly for the scalar fluctuations. We 
discuss below the modifications arising from the inflaton 
potential and the mixing with metric degrees of freedom. 
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The tensor fluctuations of the massless graviton obey ac¬ 
tually the same evolution equation as for a massless scalar, 
up to an overall normalization of the inhomogeneous term. 
Eqs. (IMl) . (11351) . (11381) can therefore be employed for the 
tensor power spectrum. Up to overall normalization eq. 
([ 2 ]) describes then the graviton spectrum during inflation. 
This tensor power spectrum remains unchanged after the 
end of the inflationary epoch as long as the wavelength of 
the mode stays far outside the horizon. Once a tensor mode 
enters again the horizon in the radiation or matter domi¬ 
nated epoch the solution of the mode function will change 
qualitatively, leading to a fast decrease of the amplitude. 

3. Inflaton potential 

Realistic models of inflation are often described in terms 
of an effective potential V((p) for the inflaton. This gener¬ 
alizes the effective action for the scalar field to 

T = J ei^^d^ipd^ip + V{ip)^ . (166) 

The second functional derivative is extended to (e = i) 

r(2) = ia^D + ia^rr?{(p), (167) 

where we have not written explicitly the unit matrix 5(x — 
y). The mass term. 


(168) 

is evaluated for the homogeneous background scalar field 
(p{r]) as determined by a solution to its field equation 


Depending on the size of the dimensionless ratios 


9 9 9 ~ 9 9 9 

2 am m m am 


(171) 


there are different regimes for the time evolution of Wk- We 
concentrate here on the case of constant (or slowly varying) 
ifi^. The evolution equation for Vk takes the form 


{d^-{2 + 1 ^-rh^)n^+ k^)v = 0. (172) 

We therefore, replace in the last term in eq. ()138|) the factor 
2 + 1 /hy 2 + 1 / — and find for the leading behavior for 
y -t 0 


lim 6 (j/) = 6 o?/ 1 + 2 , 
y^O 

with 


2i/ ^ /i 4m2 \ 

= + " (3 +12)2 J • 


(173) 


(174) 


In consequence, the equal time correlation function obeys 


k^G{k,y) 


\b^?H^ ( k \ 

2 \aH) 


(175) 


For constant v and w? this result is exact. It holds for arbi¬ 
trary values of v and . For cr > 0 the Fourier transform 
of G'(fc, rj) exists. We obtain again 

k^G{k,a —>■ oo) = = 1), (176) 

with |5op and G adapted to eg. (11721) . The large time be¬ 
havior of the correlation function for a massive scalar field 
is further discussed in appendix D. 


a +—{(fi) = 0. (169) 

dtp 

We emphasize that eq. dUZl) holds for an arbitrarily 
non-linear inflaton potential V{p), thus including interac¬ 
tions. It is by no means restricted to a quadratic potential. 
The mass term m? may depend on time through the time 
dependence of the background field. For ^ 0 the matrix 
is invertible if we restrict the function space according 
to our conjecture on the behavior of the analytic continua¬ 
tion or by imposing initial conditions. A constant field (/Jq 
no longer obeys T^^^yjo = 0. A mass term is a simple way 
to regularize r*^ 2 )^ 

The time evolution equation for the correlation function 
in momentum space remains valid, now with 

{d^+ 2TLdrf + k'^ + a'^m'^)wk{ri) = 0. (170) 

Let us consider models where for 77 —>• — oo, a —>■ 0 the mass 
term remains bounded or increases less fast than a~^. It 
can then be neglected for 77 —> — 00, such that the nor¬ 
malization of G can be computed from the “initial value 
condition” G{k,T],r]') —>■ Go{k,r],ri') for 77 , 77 ' —>■ —00. As 
time increases, however, a^m? increases and will finally 
overwhelm For a suitable potential inflation ends and 
the geometry extends to the range —00 < 77 < 00. 


4. Primordial cosmic fluctuations 

At this point we can make contact with the usual de¬ 
scription of primordial scalar fluctuations in inflation. In 
the presence of gravitational fluctuations one is interested 
in the fluctuations of the quantity 

TZk = -—Wk, (177) 

+ 

evaluated in zero spatial curvature gauge. This quantity 
is simply related to gauge invariant scalar fluctuations and 
has the property that it remains constant once a given k- 
mode has left the horizon. We can therefore employ the 
correlation function for 72., 

cH^ 

k^Gn{k,y^0) = ^\y = l, (178) 


with c a factor close to one which is the analogue of 
in eas. (11621) . (11761) The scalar spectral index Ug is defined 
by 


Us - I 


g]^ln [k^Gnik,a^oo)] 
^ dlnH 


Ai/ 2p 

1 + 1/ (1 -I- i/)H(p 


(179) 
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All quantities in the last two lines of eq. (I179|) are evaluated 
for 2 / = 1. 

The fluctuation amplitude is directly observable in the 
CMB. For the power spectrum 

P{k) = Gnik,r]^0) (fSO) 


one uses conveniently the normalization 

A^{k) = ^Pik) = Myr^-\ (i8i) 

The values observed by the Planck satellite are 

= 2.2 • 10"®, ks = O.ObMpc-^. (182) 
Comparing with ea. (ll78l) yields, approximating c = 1, 

= ^\y = r (183) 


with ps = ks/aH. 

For small iz and we recover the standard results of 
slow roll inflation. For a sufficiently flat potential eq. (I169p 
can be approximated by 


dp dV 


(184) 


and the gravitational field equation for Einstein gravity 
becomes (with M the reduced Planck mass) 

3M^H^ = V{(p). (185) 


This allows us to express v and fh? in terms of the standard 
slow roll parameters e and t] 


_H _ 1 dV dp _ kdlnVy _ 

772 dp dt 2 \ dp ) 

(186) 

and 


With 

0 = 

(p = 

one finds the standard result 

2 

ris = 1 + 6:z + -rr? = 1 — 6e + 2?7. (189) 

o 

Also for the amplitude the slow roll approximation yields 
the usual result 


3M2 52 F ^ 

-3 = 377. 


V dp^ 


(187) 


1 dV 
3 H dp ’ 
I'dV H 
3 dp 


m^p 


(188) 
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We conclude that the solution of the propagator equa¬ 
tion according to our conjecture on the behavior of the 


analytic continuation yields a scaling correlation that re¬ 
produces the usual discussion of fluctuations in inflationary 
cosmology. Results are often even exact and not only ap¬ 
proximative. Nevertheless, we recall that the realization 
of the scaling correlation is not guaranteed. It has to be 
demonstrated that is is reached before horizon crossing, ei¬ 
ther by virtue of a non-linear equilibration process, or by 
restrictions on initial conditions. If not, the power spec¬ 
trum will be modified as we discuss next. 


VIII. Impact of initial conditions on pri¬ 
mordial fluctuations 

In this section we turn back to the initial value prob¬ 
lem and ask more quantitatively: how important is the in¬ 
fluence of initial conditions for the primordial fluctuation 
spectrum in inflationary cosmology? We work here within 
an approximation of the effective action which keeps for 
the inflaton a covariant kinetic term and a potential. Pos¬ 
sible non-linear symmetrization and dissipation effects, as 
discussed at the end of sect. El are omitted. We investi¬ 
gate rather general initial conditions, that can be set many 
e-foldings before the observable fluctuations cross the hori¬ 
zon. Within this setting we find a substantial possible de¬ 
pendence of the fluctuation spectrum on the initial condi¬ 
tions. If de Sitter symmetry is not imposed on the initial 
conditions rather arbitrary fluctuation spectra can be ob¬ 
tained for suitably chosen initial values. The inflaton po¬ 
tential no longer dictates the shape and amplitude of the 
primordial fluctuation spectrum. It rather distorts slightly 
a preexisting initial spectrum. This underlines that predic¬ 
tive power for inflationary cosmology either needs strong 
enough non-linear symmetrization and dissipation effects, 
or a detailed understanding of the selection of initial con¬ 
ditions. 

Our conclusions on the strong impact of initial condi¬ 
tions seems at first sight to differ from the suggestions in 
ref. 0 that de Sitter geometry leads to an attraction to¬ 
wards universal physics even in the absence of non-linear 
terms. While the underlying equations are similar, the dif¬ 
ference in conclusion originates from three observations. 
For a massive scalar field, we find indeed for every /c-mode 
an asymptotic attraction to a universal form, in accordance 
with ref. 0 . This asymptotic solution results from the ex¬ 
ponential decay of one of the mode functions relative to 
the other. It does not fix the freedom of an independent 
amplitude for each fc-mode. Still, the correlation function 
becomes unique asymptotically if one imposes scale or de 
Sitter symmetry. For realistic scalar masses in inflationary 
cosmology, however, this asymptotic attraction sets in only 
a long time after the fluctuations have crossed the horizon. 
At the time of horizon crossing, which is relevant for the 
observable fluctuation spectrum, the information from ini¬ 
tial values is still present - memory is not yet lost. 

The second and third points concern the selection of ini¬ 
tial values for a massless scalar field. Only a small subset of 
possible initial conditions is investigated in ref. i- They 
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are based on a particular assumption motivated by infrared 
finiteness. We do not impose such restrictions since the in¬ 
frared behavior is not a severe problem for physical observ¬ 
ables. Formally, it can always be cured by changing the ini¬ 
tial spectrum for extremely small k that are far outside the 
observable range. Since in our approximation the different 
fc-modes follow separate evolution equations, a change of 
the spectrum in the extreme infrared does not affect the 
investigation for the observable fluctuations. As a third 
difference ref.Q considers pure quantum states, while our 
investigation admits mixed states as well. The setting in 
ref. @ is similar to ref. 0 . Also ref. Q suggests a partial 
insensitivity of observables from the choice of initial con¬ 
ditions. This is again due to the selected observables and 
restricted initial conditions. 

1. Time evolution of correlation functions in 
inflationary cosmology 

For a scalar field in a given geometry and with arbitrary 
inflaton potential one finds for the functions A, B and C in 
the general form of the second functional derivative of the 
effective action cm 

A=l , B = P+ a^m^ , C = -H. (191) 


a small effect. For realistic inflationary models one has 
w? ^ 1 before horizon crossing. With m} jiA of order 
one only a long time after horizon crossing the role of the 
mass is rather secondary for the issue if initial information 
is lost. 

For most of the time before horizon crossing \h/u\ is also 
small. In this range one has the general solution 

Gipip = a{k) + P{k)s\n{2u +5{k)), 

Gttip = /3(fc) cos (2u + i5(fc)), 

= a(k) — P{k) sin (^2u + S{k)^. (196) 

The free functions a{k),l3{k) and 6{k) specify the initial 
conditions. If initial conditions are set early enough, a{k) 
coincides with a{k) in the general solution (IMl) . Similarly, 
/3(fc) and S{k) are related to /3(fc) and 7 (fc) by 

/3(fc) cos(2m) — 7 (fc) sin(2u) = /3(fc) sin (2u -I- 6{k)). (197) 

It is obvious that in the range |m| ^ 1 memory is not lost. 
The scaling correlation corresponds to /3 = 0 and a = 1. 
Other initial conditions have rather arbitrary real functions 
a{k),j3{k) and 5(k), obeying only mild restrictions. 


Here irk' and % are functions of rj according to the par¬ 
ticular geometry and cosmological solution for the inflaton 
mean value <^(??). It is straightforward to solve the evo¬ 
lution equations for the correlation functions numerically. 
We employ dimensionless quantities u = krj and 


2 a^ 

Gipip = 2a^kGioio , Gttm = 2a^GTruy , Gt^tt = —G 




TTt/? ? 'TTTT 


In these units the evolution equations become 


(192) 
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^TTTT \ ~r O I 'GJ/ninn, 

V mV 

—Ctt-tt — 2 1 H— G-r^ip, (193) 
u \ ) 


with 


h = —Hr] = —T] 


dine 

dr] 


« 2 2 2 2 
, m = a T] m . 


(194) 


For de Sitter space one has 


h = 1 , = — 7T = rh^. 


(195) 


Eq. (11931) is a system of linear differential equations. In 
the limit of constant h and rh? it is independent of k. In 
this case the dependence of the fluctuation spectrum on k is 
completely described by the relations between u and ry and 
by the relation (11921) between G and G. If h and de¬ 
pend on 7] this induces an additional fc-dependence through 
h(u/k),rh'^(u/k). Typically, this type of dependence is only 


2. Initial conditions 

If we set initial conditions sufficiently far in the past, 
i.e. for |u| 1, we can use eq. (11961) for their char¬ 

acterization. Compatibility with the form of the leading 
expression Go(r, rj, r]') for the propagator in position space 
for r —>■ 0, ?7 —>■ ?]' (cf. eq. (lA.ll) ! is achieved for 

lim d{k) = 1, (198) 

k—^oo 

provided the Fourier transform of the contributions from 
d{k) — 1 and from the term ^ /3(fc) remain sufficiently 
subleading. The asymptotic form Gq can be associated 
with the Hadamard condition for a well posed boundary 
problem in position space. (In momentum space the initial 
value problem es well posed anyhow.) For d{k) — 1 ~ k~k 
the difference between G{r,r]) and Go(r, ry) is proportional 
rA“2 fQj. ^ ji; jg subleading for fi > 0. Similarly, 

for i3{k) ^ and constant <5 the contribution to G(r) — 
Go{r) remains subleading for /2 > 0. We remark that these 
restrictions concern only the limiting behavior for fc —>■ oo, 
leaving much freedom in the range of finite k. For fixed ?y 
the limit fc —>■ oo corresponds to it —>■ —oo such that the 
approximation (I196p is indeed appropriate. 

A similar discussion can be done for the infrared behav¬ 
ior. For fixed ry the limit fc —>■ 0 corresponds to u —^ 0. Now 
the terms ^ h and rh? matter for the issue of infrared finite¬ 
ness. Specifying initial conditions for some large negative 
r]in the modes relevant for the shape of the spectrum in the 
far infrared are already out of the horizon at this time. For 
these modes we can require suitable boundary conditions 
in order to have a well defined limit both in momentum 
and position space. The dynamics of such modes far out¬ 
side the horizon is, however, of no concern to the fate of 
the modes with wavelength that become observable. The 
observable modes are all well within the horizon at lym, 
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with large \uin\- Their horizon crossing occurs much later, 
when u reaches —1. In a sense, the extreme infrared modes 
only form a background for the observable modes. Their 
influence can be incorporated effectively into the effective 
action T that is relevant for the observable modes. 

Besides the conditions on the ultraviolet and infrared 
behavior which are not very restrictive for the observable 
modes, we may also require that the power spectrum should 
be positive at all times. The condition G^p^p > 0 is obeyed 
for \i3{k)\ < a{k),a{k) > 0. We will actually use the 
stronger “fluctuation constraint” (I132p which amounts to 

a{k) > 1 , |^(A;)p < a'^{k) - 1. (199) 

3. Processing of initial spectrum by geometry and 
inflaton potential 

At a time close to horizon crossing, but still sufficiently 
before, say u = —10, the initial information is still fully 
preserved. The spectrum of fluctuations at this time 
can be associated with the initial spectrum. As a given 
mode goes through the horizon, say during the interval 
—10 < M < —0.1, the details of the geometry and inflaton 
potential lead to a distortion of the initial spectrum. This 
typically violates scale symmetry since different k modes 
reach the horizon for different values of r]. Nevertheless, 
this “processing” of the initial spectrum near horizon cross¬ 
ing is a small effect since scale symmetry is approximately 
realized during inflation. Already a modest deviation of 
the initial spectrum from the “universal form” a = 1, /3 = 0 
can be comparable to the effect of this processing. In short, 
the physics relevant at horizon crossing does not “produce” 
the primordial fluctuation spectrum. It rather processes an 
“initial spectrum” that keeps memory of the “beginning of 
inflation”. The properties of the detailed geometry at the 
time of horizon crossing and the shape of the inflaton po¬ 
tential do only determine this processing. 

The situation is particularly simple for /3(fc) = 0. In this 
case every fc-mode follows precisely the same evolution as 
for the scaling correlation. Only the amplitude is multi¬ 
plied by a{k). As a result, the correlation function reads 

G{f],k)=a{k)G{r],k), (200) 

where G{rj, k) denotes the scaling correlation (1901) . In con¬ 
sequence, the amplitude of the cosmic fluctuations is mul¬ 
tiplied by a(k). This can be a factor much larger than one. 
For the spectral index one finds 


Us + Up, 

(201) 

d\na{k) 

dink 

(202) 


and Tis the spectral index for a = I. 

4. Example for initial conditions 

The initial conditions at the beginning of inflation may 
depend strongly on the particular model for the pre¬ 
inflationary era for our universe. It could be formation 


of a bubble in chaotic inflation, inflation could last forever 
since the infinite past, there may be a bounce and many 
more possibilities. Nevertheless, it seems useful to consider 
a definite example for an initial correlation function, such 
that issues of scales become apparent. 

For the beginning of inflation we consider a correla¬ 
tion function that approaches the propagator in Minkowski 
space for wavelengths much smaller than the Planck length, 
and has large classical fluctuations on length scales larger 
than the Planck length. Inbetween, there is a smooth 
crossover, with width B in logarithmic fc-space. Such a 
spectrum may be modeled by 

A 2 k 

a(fc) = I-|-— (1- arctg x)^ x = i?“^ln( —), (203) 

2 TT Kq 

and /3(fc) = '){k) = 0. For the classical fluctuations on 
large length scales we assume an amplitude much larger 
than the one for a pure quantum state, A ^ I. Thus a{k) 
varies between one at short distances and A -|- I at large 
distances. 

We need the value of x that corresponds to the ob¬ 
servable fluctuations that leave the horizon around 60 e- 
foldings before the end of inflation. This value will depend 
on the number Nin of e-foldings between the beginning 
of inflation and horizon crossing of the observable fluctua¬ 
tions. For the beginning of inflation we assume k/k^ = p, 
where p is the physical momentum of a given fluctuation 
in units of the Planck mass M. We denote by am the 
scale factor at the beginning of inflation, and ahc the one 
at horizon crossing. At the beginning of inflation one has 

— = pM, (204) 

^in 


while at horizon crossing one finds the physical momentum 



^hc 


This yields 


H^ahc 

Mam 


= e^'' 


E2. 

M’ 


and therefore 


X = B 


-1 





(205) 


(206) 


(207) 


In terms of am the crossover scale fco is given by 


k[) — amAI, 


(208) 


With eq. (12071) for x we can compute the amplitude and 
spectral index for the observable fluctuations. Assuming 
here for simplicity that the slow roll approximation is valid 
at horizon crossing one finds for the power spectrum 


A\k) 


a{k)V 
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(209) 
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There is a substantial enhancement of the amplitude for 
X A, while for x ^ A the Bunch-Davies spectrum be¬ 
comes a good approximation. For the spectral index one 
finds 


ns = 1 -I- rip — 6e -I- 2r], (210) 


with 


91nQ;(fc) 


dink 


B{1 + x‘^){— — arctg x) 


1 -1 


For a; ^ 1 this yields 


1 


=-is = 


-1 


( 211 ) 


( 212 ) 


with larger |np| for smaller values of x. Small values 
|np| <C 6e — 2ri require a large number of e-foldings Nin- 
Typically, a correction of less than 10% of the standard re¬ 
sult for ns — 1 needs for our example a rather large number 
of e-foldings between the beginning of inflation and horizon 
crossing ^ 250. This demonstrates the important role 
of an understanding of initial conditions for the quantita¬ 
tive interpretation of CMB-data in terms of the inflaton 
potential or similar quantities. 

On the other hand, cosmologies with a very long, perhaps 
eternal, duration of the inflationary epoch of our universe 
can provide for an effective approach to the scaling cor¬ 
relation in a rather simple way. For —>■ oo the effect 

of the initial conditions becomes negligible in our exam¬ 
ple. In this limit the relevant part of the initial spectrum 
is always the extreme ultraviolet tail. If this is given by 
the universal propagator for Minkowski space one finds at 
horizon crossing of the observable fluctuations the scaling 
correlation. This is in line with the general findings of 
sect. El Initial values given by the scaling correlation lead 
to a scaling correlation for all later times. For observations 
it is sufficient that this holds for the relevant range in k. 
For Nin 00 this relevant range is given by the ultravi¬ 
olet limit k ^ oo, for which the universal propagator in 
Minkowski space can be motivated by various arguments. 
We may be able to look back to the beginning of inflation. 
What we may see, under certain circumstances, are simply 
the Lorentz-invariant fluctuations in flat space. 


IX. Conclusions 

Do the observable primordial cosmic fluctuations keep 
memory of the “beginning of the universe” or the “begin¬ 
ning of inflation” ? The answer to this question is not yet 
entirely settled at present. This paper develops the con¬ 
ceptual setting for investigations that can clarify the sit¬ 
uation. It addresses questions about the “correct choice” 
of the Green’s function for de Sitter space, the “correct 
choice” of quantum vacua for time evolving situations or 
more generally the determination of a universal scaling cor¬ 
relation. 

The main unknown for a complete answer is the “equi¬ 
libration time” it takes for an effective loss of memory of 


initial conditions for fluctuations in the observable range. 
At the present stage it is even not settled if equilibration 
takes place at all. Since interactions are small during in¬ 
flation one expects, at least, that a possible equilibration 
time is very long, covering many e-foldings. Looking back 
towards the beginning of the universe, observations can de¬ 
tect signals from early stages of inflation, dating back for 
a time at least as long as the equilibration time. 

The key ingredient are time evolution equations for the 
correlation functions that permit to follow their fate for 
rather general “initial conditions”. If memory of the de¬ 
tailed initial conditions is lost rapidly enough, a unique 
propagator is selected as an asymptotic attractor solution. 
Only in this case the choice of initial conditions is unim¬ 
portant for an understanding of the cosmic fluctuations. 
Their properties will then reflect the physical situation in 
the epoch when the observable fluctuations leave the hori¬ 
zon. Observation of structure in the universe can therefore 
determine properties of the inflaton potential or similar 
quantities. If, in contrast, information about initial con¬ 
ditions still plays a role at horizon crossing, observations 
can only constrain some combination of physics at horizon 
crossing and “initial conditions”. 

The investigation of the present paper has been limited 
to a scalar field in a given time-evolving background ge¬ 
ometry. Interactions are included in form of the inflaton 
potential. We have made an approximation where non¬ 
linear effects due to backreaction and an explicit depen¬ 
dence of the effective action on the correlation function are 
neglected. In this approximation we find that the infor¬ 
mation about initial conditions is not lost - equilibration 
time is infinite. Typically, infinitely many conserved quan¬ 
tities are an obstacle to equilibration. If our approxima¬ 
tions are valid since the beginning of inflation, the observed 
anisotropies in the CMB reflect the initial fluctuation spec¬ 
trum at the beginning of inflation, with only a modest pro¬ 
cessing at the time of horizon crossing. 

For a very long duration of inflation, perhaps even since 
the infinite past, we find it likely that our approximations 
do not remain valid and some equilibration occurs, at least 
for the range of wavelengths that are observable and which 
have been in the extreme ultraviolet in the remote past. 
For this case we formulate a conjecture that a unique cor¬ 
relation function is selected as an attractor solution. This 
is similar to the approach to thermal equilibrium. We give 
a criterion, based on the behavior of the analytically con¬ 
tinued Green’s function in the infinite past, which selects 
this universal scaling correlation uniquely for a given form 
of the effective action. For free scalar fields in de Sitter 
space this singles out the Bunch-Davies vacuum, while it 
can be extended to a large variety of situations, including 
interacting theories. 

A possible equilibration needs non-linear effects beyond 
our approximation, as backreaction or the explicit depen¬ 
dence of the effective action on G. Even without a de¬ 
tailed investigation of these effects some general features 
of the approach towards a unique correlation function be¬ 
come visible by our investigation. If the correlation func¬ 
tion has reached a state with approximate de Sitter symme- 
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try or scale symmetry already before horizon crossing, the 
shape of the power spectrum is close to the one assumed 
in the standard scenarios of inflation. This concerns, in 
particular, the spectral index. Symmetry is, however, not 
strong enough to fix the relation between the fluctuation 
amplitude and the Hubble parameter at the time of horizon 
crossing. For the latter some type of fluctuation-dissipation 
relation is needed. In the presence of de Sitter symme¬ 
try the possible deviations from the standard setting affect 
only the extracted value of the Hubble parameter at hori¬ 
zon crossing. We conclude that a rather rough approach 
of the correlation functions towards an equilibrated state 
is sufficient to realize the overall picture of fluctuations in 
inflationary cosmology. This is analogous to prethermaliza- 
tion [ 5 ^ where only rough properties of correlation func¬ 
tions can reproduce some characteristic features of thermal 
equilibrium states. 

When it comes to precision and the determination of the 
inflaton potential by observation of the cosmic microwave 
background, the situation is less clear. The deviation of the 
spectral index Ug from one is already a small quantity, such 
that rather modest memory of the initial state could dis¬ 
tort it. For Us — 1 one has to compare the scaling violation 
arising from the time evolution of the Hubble parameter 
and the effective inflaton mass term (the parameters e and 
rj in the slow roll approximation) with the size of effects in¬ 
duced by memory of the initial state. Only if the latter are 
negligible the usual determination of the inflaton potential 
by measurement of Ug and the tensor to scalar ratio r is 
justified. 

Even if an approach towards a unique correlation func¬ 
tion takes place (what we believe), the issue of the time 
scale for such an equilibration process remains important. 
For a scalar theory in a given background this “equilibra¬ 
tion time” would diverge in the limit where the interactions 
go to zero. In the presence of metric fluctuations one has to 
asses the quantitative role of the gravitational interactions 
that are always present. For the interpretation of observa¬ 
tions the time scale for the approach to a symmetric state 
(“symmetrization” or “generalized isotropization” H) is 
particularly relevant. 

Symmetrization and more generally equilibration is not 
expected to affect all wavelengths equally. The extreme 
infrared scales A: —>■ 0 are already outside the horizon very 
early in inflation and therefore not subject to equilibra¬ 
tion processes. On the other hand, for the ultraviolet limit 
of the spectrum, /c —>■ 00 , one may expect equilibration 
on comparatively short time scales. If couplings are taken 
dimensionless and k/a exceeds all mass scales in the the¬ 
ory, one expects for the equilibration time a relation of 
the type teq = <j/k, with a a (possibly very large) dimen¬ 
sionless factor depending on couplings. This follows from 
dimensional analysis. (For vanishing interactions or in the 
presence of conserved quantities that obstruct equilibration 
a diverges.) In the extreme ultraviolet limit the evolution 
equation for the correlation function approaches the one for 
a free massless field in Minkowski space. It is reasonable to 
assume that for fc —>■ 00 the asymptotic solution is the free 


propagator in flat space. The factor a can be determined 
by a computation in flat space. 

If the number of e-foldings Nm between the beginning of 
inflation and horizon crossing of the observable fluctuations 
tends to infinity, the above considerations select for the ob¬ 
servable range of modes the scaling correlation according to 
our conjecture. This holds for any finite cr. For N^n —t 00 
the observable fluctuations correspond to /c —>■ 00 . Also the 
time towards the beginning of inflation diverges, such that 
equilibration becomes arbitrarily accurate. On the other 
hand, for small Nin the time since the beginning of infla¬ 
tion will not suffice to produce a strong equilibration. In 
this case the state of the universe at the beginning of infla¬ 
tion becomes observable. Inbetween, there will be a value 
such that for ^ the memory of initial 

conditions is effectively lost, while for ob¬ 

servations are sensitive to the state of the universe at the 
beginning of inflation. The main open point for an answer 
to our question concerns the quantitative determination 
of This may depend on the particular inflationary 

model. 

We have addressed the issue of the influence of initial 
conditions within the framework of the quantum effective 
action F. This offers the advantage that all effects of fluc¬ 
tuations are already incorporated in the effective action. 
The second functional derivative of F yields an exact rela¬ 
tion for the propagator and determines its time evolution. 
The complicated discussion of quantum vacua in a time- 
evolving situation and the “correct choice” of such a vac¬ 
uum are avoided in our setting. Once the effective action 
is computed or assumed, the determination of the Green’s 
function becomes an issue of “classical” field theory. Nec¬ 
essary approximations concern only the precise form of the 
effective action. 

Powerful functional methods for the computation of F 
are available, as functional renormalization or the non¬ 
equilibrium two-particle irreducible effective action. We 
hope that the relevance of initial conditions for the cosmic 
fluctuations can be settled quantitatively by such methods. 

Appendix A: Correlation function in 
position space 

In this appendix we discuss the propagator in position 
space. We first discuss the general short distance behavior. 
It is the same as for flat space. This limiting behavior for 
vanishing distance is often associated with the Hadamard 
property. We will see that the short distance behavior of 
the propagator does not fix it uniquely. All correlation 
functions discussed in this note actually share the property 
of a limiting Minkowski behavior, but differ from each other 
at larger distances. Next, we concentrate on the propaga¬ 
tor of a massless scalar field in de Sitter space, discussing 
first the general short distance behavior and subsequently 
the Fourier transform of the scaling correlation in momen¬ 
tum space. The latter does not exist due to an infrared 
divergence of the momentum integral. We also investigate 
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a proposed propagator that is well defined in position space 
and consistent with scale or de Sitter symmetry. It does 
not correspond to the scaling correlation. Even more, the 
associated initial conditions do not correspond to a mixed 
quantum state. 


1. Short distance behavior of propagator 

The short-distance behavior of the propagator in po¬ 
sition space is closely related to the one in Minkowski 
space. Indeed, the leading behavior of the Green’s func¬ 
tion G(r, rj, T]') for r —?> 0 and rj ^ rj' is given by 


47 r 2 a(r?)a(r?')(r^ — (1 — 2 ie )(77 — 77')^) ’ 

(A. 1 ) 

which replaces in the free Green’s function ( 1551 ) the factor 
by 0(77)0(77'). With 

G = Go + < 5 iG, 

= (1 + 276 ) 52 -A, 

D = DQ + 2 {l + 2 ie)ndr^, (A. 2 ) 

we observe 

i{l-it)Do[a{vi)a(jq')GQ\ = 8{jq - T^') 5 ^{f). (A. 3 ) 

This yields 

D 5 iG = ^^Z?o[a(??)Go] -f?Go 

= (l + 2 ze)('H 2 + 5 ^'H)Go. (A. 4 ) 

Thus 5 Gi is determined by an inhomogeneous differential 
equation in terms of Gq. The r.h.s. is less singular than 
the ^-distribution in eg. ( 1301 ) . The leading singularity is 
therefore taken care of by Gq . 

For the particular geometry of a radiation dominated 
universe (77 > 0), 

TL = — , a = Cl77 = C2\/t, (A. 5 ) 

V 

the r.h.s. of eq. (IA. 4 I) vanishes and Go is the exact propaga¬ 
tor. For the general case we may explore near r = 0 ,77 = 77' 
a propagator of the form Gq + Gi, with 

^ / r + (77-77') y 

^ 2 >2 tt‘^ a{r])a{ri')r — {rj — rj')) ’ 

K{r],ri') = "^(77) — "^(77') + J dfj'H^if}). (A.6) 

rj' 


We note that Gi is less singular for r —^ 0 ,77 —^ 77' than Gq . 
(We have set e = 0 for Gi, and similarly on the r.h.s. of 
eq. (lA. 411 .1 Inserting eq. (|A.6|) into eq. (IA. 4 I) yields (r 0 
or 77 77') 


S,G = Gi + <52G, 

D{ 62 G) = {n^ + dr,n)Go - DGi 

= {{n^ + d^n)K - 2ndr,n - d^n) 


(A. 7 ) 


lnfl±2^y 

\r-v+ri J 


S2Tr^a{ri)a{ri')r 


where H stands for "^(77). As compared to D{SiG) in eq. 
(IA. 4 I) we see how the degree of divergence of D{S2G) for 
r —>• 0 ,77 — ^ 77' has decreased further. We observe that Gi 
vanishes for rj = ij'. The form of Gi is not unique, however 
- the contribution 62G may have a similar divergence for 
r —>■ 0 , rj ^ rj'. Furthermore, DGi contributes for r = 
0 ,77 = 77' to an inhomogeneous term. This contribution has 
to be canceled by 52G. 

A unique short distance behavior in position space, as 
given by Gq, does not imply that the limit /c —^ 00 in 
momentum space is unique. We will encounter examples 
where G{k —>■ 00,77,77') differs from [20(77)0(77')/;:]“^, while 
in position space Go is approached for r —>■ 0, 77 — >■ 77'. 


2. Short distance propagator in de Sitter space 

Let us consider the particular geometry of de Sitter 
space, 

•H = -- , o = -7^, (A.8) 

77 Hor] 

where 


1^(77,77') = 2 (H( 77 )-H( 77 ')). (A. 9 ) 

(The range of 77 is here —00 < 77 < 0 , with t = 
— ln(—iJo 77 )/iLo-) Eq. (IA. 7 I) simplifies to 


D{ 62 G) = - 


STT^r 


In 


r + (77 — 77') 
r — {jj — 77') 


For this case we may investigate 

5 ^G = G2 + 63G, 

with 


G2(r, 77, 77') = - 


nivmvl 


1677^0(77)0(77') 


In 


(A.IO) 


(A. 11 ) 


[r^ — (77 — 77')2]2 


+ !^ln 


(r + 77 - 77')^ 
(r — 77 + rj'y 


(A. 12 ) 


For the particular case of de Sitter space one finds 

DiSsG) = 0. 

The combination G = G0+G1+G2 is an exact solution of 
the propagator equation (I 29 p , including the inhomogeneous 
term. The propagator 

G(r, 77 , 77 ')= (A. 13 ) 

Hq f W 

4jr2 ^2 _ 


T9 A — In 
)2 4 


[y.2 _ (y, _ y,7)2]2 


involves an undetermined scale tq which corresponds to 
the possibility of a constant shift in G. The equal time 
propagator for 77' = 77 takes the simple form 


G(r,77) = 



(A. 14 ) 


The dependence on 77 — 77' is independent of the scale tq. 
In particular, one has 


G2(r, 77, 77') - G2(7', 77, 77) 


IGtt^ 

77-77' 


In 


(A. 15 ) 

(r + 77 - 77')^ \ ) 
(r- 77 + 77')2y / ■ 
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The function G{r,T],r]') in eq. (IA.13I) solves equation 
(Ei) - the contributions of Gi and G 2 to the inhomoge¬ 
neous term cancel. Nevertheless, G{r,rj^f}') may not be 
considered as an acceptable propagator for de Sitter space. 
In particular, G(r, rf) turns negative for large r and the 
Fourier transform is not defined. Furthermore, G(r, 77 , 77 ') 
is not invariant under the symmetries of de Sitter space. 
One should therefore consider G(r, 77 , 77') only as a good ap¬ 
proximation for small r and (77 — 77 '). We will turn to this 
issue in sect, appendix E. 


This coincides with G2{r, rj, 77 ') — G 2 (r, 77 , 77 ) in eq. (IA.15I) . 
For the equal time part of G 2 one has {x = kr) 


G 2 (r, 77 ) = G 2 (r, 77 , 77 ) = - 


Hl 

47r^e^r J k‘^ 

0 


00 

Hq f dx 


dTT^e^ J x"^ 

0 


sin a;. 


sin(fcr) 

(A.20) 


3. Fourier transform of de Sitter propagator to 
position space 

We next turn to the conjectured scaling correlation. The 
de Sitter propagator (TOj) does not admit a Fourier trans¬ 
form to position space. If we look at the different pieces in 
eq. (with / = Hq ) we see that the symmetric part of 
G 2 diverges for fc —^ 0 as k~^ such that the Fourier integral 
(U) becomes infrared divergent. Let us consider the form 
of the different pieces Gq, Gi and G 2 in position space sep¬ 
arately. The Fourier transform of Gq reads for arbitrary 
e = cos 175 - 1-7 sin 75 and rj > ij' 


Go (g ??,»?') 


I 

4 TT^a{r])a{r]')(^r^ + 6^(77 — 77 ')^) ’ 


(A.I 6 ) 


in accordance with eq. dSil). For Gi one obtains 


Gi(r,?7,77') = 
iHq 


2e 

1 


0 ( 77 ) 0 ( 77 ') 


giferg-efe(r)-r)') 




IbTT^er \a{r]) a{r]') 

r'^ + {rj — 77')^ -I- 2 r{r) — rf) sin 79 
7 "^ + (77 — 77 ')^ — 2 r(r] — rf) sin 73 


(A.I7) 


+ 


X In 
Hq 

Sir^er \a{r]) a(r]') 


arctg 


2 r (?7 — 77 ') cos 79 
(77 — rj')"^ — 


For 79 = 7 r/ 2 , e = i the first term in eq. (IA.I7|) equals Gi 
in eq. (IA. 6 [l . whereas the second term vanishes. 

We next consider the difference between G 2 (r, 77 , 77 ') and 
its equal time counterpart G 2 ( 7 ', 77 , 77 ), 


tt 2 r pikr 

G2(r, 77, 77') - G2(r, 77, 77) = ^ - !)• 

(A.I 8 ) 

This is a well defined Fourier integral for Re{e) > 0. For 
Minkowski signature e = i the real part reads 

■Re[G 2 (r, 77 , 77 ') - G 2 (r, 77 , 77 )] = 

((r--77-f 77')2) 

r -I- 77 - 77 ' / (r -I- 77 - 77 ')^ \ I 

(r - 77 -h 77')2 y J ■ 


This integral is divergent, being formally independent of r 
and 77 . Due to the lack of a well defined Fourier transform of 
G 2 the propagator in position space can only be determined 
if r( 2 ) is regulated. We will turn to this issue in appendix 
D, E. 


4. Well defined propagator in position space? 

We do not consider the absence of a well defined prop¬ 
agator in position space as a very serious problem. For 
example, it disappears if we regularize the propagator by a 
small scalar mass. Also de Sitter space is presumably not 
realized exactly in nature, and modifications of the back¬ 
ground geometry may solve the issue. Furthermore, the 
universal scaling correlation may not be realized in the far 
infrared. Nevertheless, it is reasonable to ask if a well de¬ 
fined propagator in position space can exist if we relax our 
condition based on analytic continuation and turn back to 
the more general family of propagators (15^ . We will take 
here a, /3 and 7 independent of k in order to be consistent 
with the full continuous symmetry of de Sitter space. 

For general a and /3 the problem will persist since the 
infrared behavior for fc —>■ 0 is the same as for the propa¬ 
gator O • There is, however, the special choice a = /3 
(cf. eq. (1851) 1 for which the infrared divergence is canceled. 
Indeed, we find that for the propagator 

G(fc,77,?7') = wl{g)wl{r]') (A.21) 

+ 2 (^fc id') + wf (77)77;^ (77')) 


the Fourier transform exists. For e = 7 it reads 


f I 1 

\r 2 — (77 — 77')2 r 2 — (77 -I- 77')2 

^ gg f - iv + v')^f \ 

167r2 y _ j^/^2^2^ 


(A.22) 


For the short distance behavior r —^ 0 ,77 —>■ 77 ' the 
Fourier transform of wff(rf)w^{ri') is given by G in eq. 
(IA.13I) . The short distance behavior of the transform of 
G/3 = {wff {rf)w'f {-q') + wff{q)w'l{q '))/2 obtains therefore 
as G — G, i.e. 


Hi 


° (A.23) 


r 
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where ?? = (?? + f?0/2- The role of Gp is subleading, its 
dominant effect being to remove the free scale ro from G. 

In contrast, we observe that G(fc, rj, t]') has for small k a 
behavior completely different from G{k,r],ri') in eq. (1^ . 
The equal time correlation reads 


Using dimensionless variables and constants 


9 

u 


— G i^ipk , 


= krj = — 


h = 
k 

aHo 


, Q = 


fc4 ’ 


(B.4) 


G{k,ri) = ^^(l + cos( 2 A: 77 )) + ^(l-cos( 2 fc??)) 

+ -^ sin(2fcr7), (A.24) 

with a limiting behavior for modes outside the horizon k <C 
Hoa, \kri\ < 1 , 

The Green’s function for fc —?► 0 looks completely differ¬ 
ent from the behavior G ^ i?g/(2A:^) for the de Sitter 
propagator O- On the other hand, at horizon crossing, 
kja = Ho,kri = — 1 , one has 

G{kTj = - 1 ) = ^{1 - sin( 2 )) = 0.091^. (A.26) 

This differs from the de Sitter propagator (TOl) by a sub¬ 
stantial factor around ten, but leads to a spectrum with the 
same shape. It is the amplitude of fluctuations at horizon 
crossing that matters for the observable power spectrum, 
see sect. lVIll By virtue of the symmetries of de Sitter space 
the difference between the two propagators G and G does 
not concern the shape of the spectrum. It only concerns the 
precise relation between the fluctuation amplitude and the 
Hubble parameter. We will discuss the propagator (IA.2I|) 
in more detail in appendix E. 

Also the part ~ 7 in eq. (15^ is not infrared divergent 
and its Fourier transform exists. It may be added to the 
propagator G. 


Appendix B: Solution of evolution equation 
for equal time correlation function 


In this appendix we discuss the general properties of 
the solution of the evolution equation for the equal time 
Green’s function. We concentrate on a massless scalar in 
de Sitter space. Our starting point is eq. (llOSp . 

We can use the conserved quantity ac in eq. (II08I1 in 
order to express Gt^tt in terms of G^p^p and 

G,, = (gI^ + G-i, (B.I) 

such that 

d^G^^ = {Gl^ + c^H^ii^)G-l + ‘^G.p - k^Gp,^. (B.2) 

The remaining two equations (I105|) can be combined in a 
non-linear second order differential equation for G = G^p, 


QlG - IdM - + 2fc^G - = 0 . (b.3) 


G 


G 


the evolution equation reads 


du 

dh 

du 


h, 

2h „ Qu'^ 

-2(7 -|- -—I- —— 


(B.5) 


The dependence on the integration constant Cq or Q can 
be absorbed by a rescaling 


V\Q\ 



(B. 6 ) 


In terms of g and h eq. (IB.51) keeps the same form, with Q 
replaced by 1 (or —1 for negative Ca). The numerical so¬ 
lution shown in Figs. O |4] suggests an oscillating approach 
to the scaling solution 

9 = ,h = u, (B.7) 

which corresponds to eq. (I109|) . For various different initial 
conditions for g and h at some initial Uin (taken for FigslSl 
[Has Uin = — 10 ) the memory of the initial conditions seems 
at first sight to be partially lost for |u| < 1 . 



FIG. 3: Time evolution of g as function of u, according to eq. 
diu with Q = 1. For different initial conditions at u = —10 we 
observe oscillations around the scaling solution (central curve). 
Horizon crossing corresponds to u = — 1. 


In order to investigate this issue more closely we intro¬ 
duce functions 


2g h 

9 1 I 2 ’ ^ ’ 


1 -I- 


such that (for Ca > 0 ) 
r(^-ff), 


dg 

2u 

du 

I - 1 - 

dh 



du 

u^ 


(B. 8 ) 


(B.9) 




+ 


{l + u^)g {l+u‘^)g 
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FIG. 4 : Time evolution of h as function of u. The straight line 
is the scaling solution. 


One realizes the fixed point for h = g = 1, which corre¬ 
sponds to the scaling solution. For any given mode with 
a given fixed k the infinite past corresponds to large nega¬ 
tive u. In the regime ^ 1 one finds oscillations around 
the scaling solution with a damped amplitude. Indeed, for 
^ I we can approximate 


du 


= -{h-g) 


du 


= “( T-5 ) + 


9 


(B.IO) 


During the oscillations h and g are of a similar order of 
magnitude and we employ 


h = - , 6 g = g-1, 


(B.ll) 


where 


du6g = 2h, 
dyh = 


^ -{l + 5g)+ ^ 


l + Sg 


1 + Sg 


(B.12) 


Near the fixed point at 6g = 0,h = 0 we linearize the 
second equation (IB.12|) . 

duh = -2Sg, (B.13) 

and observe the oscillation with period Aw = tt, 

6 g = Cg sin(2u + ipg) , h = Cg cos(2u -|- ipg). (B.14) 


During the oscillation epoch the size of g and h decreases 
^ reflecting the Hubble damping. Thus ^ 

\/(o?k),GTrip ^ l/o^, with much larger in amplitude 
than the value (I109p . Except for the increase of a this 
behavior resembles the one in flat Minkowski space. 

Once gets of the order one and smaller the oscillation 
stops and one encounters an overdamped behavior. Con¬ 
sider the combination 


r = 



(B.15) 


which obeys 


df ^ ^2 u* f 1 A 

= r- -b 1 

du 1 + u^ 1 + u‘‘ \g‘‘ ) 


(B.16) 


The solution for small 


f = CrU , h = (1 + CrU)g 

implies that the running of g stops for it —>■ 0, 
= 2 u{h - g) = 2crU^g, 

according to 


g = go exp 



(B.17) 


(B.18) 


(B.19) 


Numerically we observe indeed that g stops for it —>■ 0 at 
a value go which depends on the initial conditions. The 
asymptotic behavior for /c —?> 0 is given by 

Gw = ■ (B.20) 

The appearance of the factor goy/c^ demonstrates that the 
memory of initial conditions is not lost. 

The apparent approach of the numerical solution to the 
scaling solution visible in Figs.[31l3]is mainly due to a com¬ 
bination of Hubble damping with the simple fact that G^p^p 
approaches a constant for it —>■ 0. For the initial values 
shown this constant turns out to be only moderately dif¬ 
ferent from the one in the scaling solution. The detailed 
memory of the initial conditions is not lost for the general 
solution. Since the integration constants go and Ca can de¬ 
pend on k, the spectrum of fluctuations needs not to be 
^ k~^. It retains information about the initial conditions. 
The same holds for the fluctuation amplitude. 

The structure of partial fixed points and the qualita¬ 
tive behavior of deviations from them remains the same 
for more realistic inflationary geometries as the ones dis¬ 
cussed in sect. IVHI This demonstrates analytically that 
there is no loss of memory of initial conditions in our ap¬ 
proximation of a derivative expansion. 


Appendix C: Fluctuation-dissipation 
relation 

In this appendix we briefly discuss some aspects of the 
fluctuation-dissipation relation. The antisymmetric part of 
G> corresponds to the spectral function p 

p{k,r],g') =2iGa{k,r],g'). (C.I) 

For a unique Green’s function the symmetric part Gg is re¬ 
lated to p. Since Gg describes the statistical information of 
the correlation function while p contains information about 
decay rates of unstable particles, this connection is called 
a fluctuation-dissipation relation. For the most general so¬ 
lution of the propagator equation Gg contains free param¬ 
eters. Memory of these parameters has to be lost when the 
propagator approaches a universal form. The “equilibra¬ 
tion” can therefore be seen as the dynamical establishment 
of a fluctuation-dissipation relation. 
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Let us investigate this relation for the scaling correlation 
in de Sitter space where 


2 I , 1 sin(fc (?7 - 77')) 


P = Hq{ (W + 


,cos{k(jq-i)) 

-(.p-v) -p- 


r2 I / /, 1 sC0s(fc(?7-77')) 


Gs - (rj-rj' + —, 


2 k 


sin(fc (?7 - 77')) 

-2F- 


We introduce 

1 


(C.2) 


V= Tiiil + v')^ 'nr=r]-r]', r]r]'= r]'^ - (C.3) 


and make a “finite interval Fourier transform” with respect 
to the time difference rjr 

e 

p{k, u},fj)= J drjre'-‘^'^^p{k, fj, 77 ^) (C.4) 


and similarly for Gg The restrictions rj < 0, rj' < 0 for de 
Sitter space require e < —2r]. One finds 


with 


P = 


f- 2 , 1 1 52 


ik 
_ 1_5 
k dk 


V 


P 4 5P 

{Se{uj + k) - Se{uj - kj), 


" 2 k r k^ “^4 9^2 


1_5 

k dk 


{Se(u; + k) + Se{uj - k)), (C.5) 


s.i^±k) = L 

ZTT J 7T{uJ ± k) 


(C. 6 ) 


One has the normalization 


dujSe{uj ± fc) = 1 , 


(C.7) 


one can easily reconstruct Gg and p from 

G>{k,uj,fj) = ^ (^fk^ + 1 + U^-k). 

(C.IO) 

For {oj — kY » 1/e^ the function 5g{oj — k) oscillates fast, 
with period 27r/e and envelope decaying ~ l/(7r(|w — fc|). 
In this range the leading expressions for the derivatives are 

df.5e(pJ — k) = —e^Se{uj — k), 

dkSe{uj — k) = —eSg(uj — k+—'^. (C.ll) 

For fee » 1 the leading expression is 

7 r/y 2 

G> = ^{f)^-e^)Sg{uj-k), (C.12) 

while for /ce <C 1 one finds 

G^ = + l)Sg{uj-k). (C.13) 

For fixed e and k the limit 77 —>■ — 00 yields the flat space 
result 

«=■“) 

In the opposite limit (w — fc)^ <C 1/e^ one can expand 

, 5 ,( 0 ; -k) = ^(^l- ie2(cc - kf + ..)j , (C.15) 

such that 

— {2ffk^{uj‘^ — k^) + 2a;^ — fc^}. (C.16) 

For arbitrary initial conditions the fluctuation- 
dissipation relation (IC.5I) for the scaling correlation 
is not realized. Equilibration can thus be monitored by 
the degree to which the fluctuation-dissipation relation is 
dynamically realized. This corresponds to the dynamical 
realization of the particle number distribution in the 
approach to thermal equilibration. 


and for e(uj ± fc) —>■ 00 the function 6 ^ becomes the d- 
distribution. With 

G>{k,uj,fj) = Gg(k,uj,fj) - -p{k,uj,fj), 
G<ik,uj,fj) = Gg{k,uj,fj) + ^p{k,uj,fj), (C. 8 ) 

and the relations 

G<{k,uj,fi) = G>{k,-uj,fj), 

Gg(k,-uj,fj) = Gg{k,u>,fj), 

p{k,-uj,fi) = -p{k,uj,fi), (C.9) 


Appendix D: Massive scalar propagator in 
de Sitter space 

The late time asymptotic behavior of the correlation 
function differs qualitatively between a massive or a mass¬ 
less scalar field. For a massive field the fast decrease of one 
of the mode functions leads to an effective loss of memory 
once the mass term becomes important. This loss of mem¬ 
ory occurs already within our approximations, without the 
need of additional non-linear effects. We discuss this issue 
here for a de Sitter geometry. The qualitative behavior is 
similar for other homogeneous and isotropic cosmologies. 
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1. Scaling correlation 

The propagator for a massive scalar in de Sitter space 
can be obtained as a special case of the discussion in sect. 
I VIII Indeed, taking in eq. (11721) the value v = 0 and con¬ 
stant m ?, one obtains the mode equation for the correlation 
function of a free massive scalar in de Sitter space. In eq. 
(jl74l) one has 


cr = 3- \/9-4 to 2. (D.l) 

For > 0 one has cr > 0 such that the correlation function 
decreases as the scale factor a increases. Ea. (ll75|) describes 
in this case the asymptotic behavior of the propagator of 
a massive particle in de Sitter space. 

This result is consistent with a simple direct computa¬ 
tion. For k/a m we may use the field equation for a 
homogeneous scalar field perturbation in the Robertson- 
Walker metric, 


{d^ + 3Hdt + m^)w = 0. (D.2) 

For the variable s = Ht = — ln(—i/r;) = Ino it takes the 
form 


{d1 + 3ds + m^)w = 0. (D.3) 

For the general solution 

u, = -bC 2 e^"®, (D.4) 


with 


Ai /2 = ^ (-3 ± \/9 - 4to2^ 


(D.5) 


the leading behavior for s —>■ oo is given by the larger 
eigenvalue (with -|- sign in eg. (ID.SI) !. The result 


licp = |cip exp |(-\/9 — — 3)s| 




(D. 6 ) 


agrees with eqs. (11751) . (ID.II) . 

For > 0 the momentum dependence of the equal time 
scalar Green’s function in de Sitter space obeys for fc —>• 0 


G{k, -q) - (D.7) 


where we use fh? <C 1 in the second expression. This cures 
the infrared divergence in eq. (IA.20I) . The Fourier trans¬ 
form of the de Sitter propagator exists for arbitrary small 
positive . Such a mass term therefore regulates the prop¬ 
agator. 


2. Time evolution 

In contrast to the massless case the scale invariant prop¬ 
agator of a massive scalar field in a de Sitter geometry 
is attracted towards a universal propagator for asymptotic 
time [H. Out of the three integration constants a, /?, 7 only 
one survives effectively. De Sitter symmetry or scale sym¬ 
metry imply then that the spectrum is unique up to an 
overall constant. This behavior is in agreement with the 


findings of ref. 0 . Without these symmetries one still has 
one free integration constant a(k) for every fc-mode. 

The (partial) loss of memory is due to the fact that both 
eigenvalues Ai and A 2 for the general solution (|D.4p are 
real and negative. The component with the larger negative 
eigenvalue (A = —(3/2 -I- \/9 — 4m^l2)) dies out rapidly, 
such that a unique asymptotic solution for Wk remains. 
The same holds for the general solution to the time evolu¬ 
tion equation for the correlation function (11051) . We note, 
however, that this attractive behavior sets in once 
is of the order m^. For small this happens only once 
modes are already outside the horizon. Before the term 
dominates in eg. (11721) and the solutions are oscillating 
similar to the solution for a free scalar field. For realistic 
inflationary cosmologies one has <C 1 and the asymp¬ 
totic attractor property plays no role for the memory of 
initial conditions for the observable fluctuations [^ . 

We illustrate this behavior in Figs. [SUH] which show the 
evolution of G as a function of m = kr]. We plot two ini¬ 
tial conditions both for a massless (blue and green) and a 
massive (red and orange) scalar field. One of the initial 
conditions corresponds to the scaling correlation (blue and 
orange), the other to a more arbitrary initial state (red 
and green). For u < —5 the effect of the mass term is 
completely insignificant, as visible in Fig IS] that shows G 
in units of k~^a~'^. In FiglSIwe focus on the behavior after 
horizon crossing {u ^ —1) and plot G in units of HQ/k^. 
At horizon crossing, u = —1, the effect of the mass is vis¬ 
ible for our choice of = 0.5. The asymptotic behavior 
for u —^ 0 shows that for the massive scalar the correlation 
function becomes independent of initial conditions, in con¬ 
trast to the massless case. For purposes of illustration we 
have taken a rather large mass term, fh? = 0.5. For the 
smaller masses that are characteristic for realistic inflation 
both the separation between massive and massless propa¬ 
gators, and the approach of the massive propagator to the 
universal asymptotic form, set in for u much closer to zero. 
The difference between the massive and massless scalar at 
horizon crossing would not be visible in the plot. 



FIG. 5 : Time evolution of correlation function G in units of 
as a funtion of m = kq. The initial conditions corre¬ 
spond to the scaling correlation (smooth curves) and another 
(rather arbitrary) initial state (oscillating curves). The differ¬ 
ence between the massive and massless case are only barely 
visible, except for the final increase for u —>■ — 1. 
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FIG. 6: Time evolution of G in units of H^/k^, for modes out¬ 
side the horizon {u ^ — 1 ). Initial conditions are the same as 
for Fig. [ 5 ] The two curves for the massive scalar field approach 
the same asymptotic propagator, independent of initial condi¬ 
tions. In these units this propagator vanishes for u —>■ 0 . For a 
massless scalar field the propagator approaches a constant that 
depends on the initial condition. 


Appendix E: Symmetries 


In this appendix we discuss the consequences of symme¬ 
tries. For a de Sitter background the scalar correlations 
that exhibit the symmetries of de Sitter space are partial 
fixed points of the time evolution. A possible “symmetriza- 
tion” would describe the approach to these partial fixed 
points. We are therefore particularly interested in correla¬ 
tion functions with de Sitter symmetry. 

General covariance dictates the transformation proper¬ 
ties of the correlation function under coordinate trans¬ 
formations. This restricts the functional dependence of 
G(x, y) on the coordinates x and y. In case of isometries of 
the background geometry one may require that the prop¬ 
agator is invariant under the symmetry transformations. 
Isometries can then impose rather severe restrictions on G. 
While Poincare symmetry fixes G uniquely for Minkowski 
space, this is not the case for the symmetries of de Sitter 
space. 

1. Covariance under general coordinate 
transformations 

The propagator in position space should not depend on 
the chosen coordinate system. A general coordinate trans¬ 
formation changes both the vierbein e^(x) (and the associ¬ 
ated metric g^u) and the scalar field g}{x). An infinitesimal 
transformation of the scalar translates directly to the cor¬ 
relation function which transforms as a scalar bilinear 


5(p 


6 Gix,y) 


-^'"{x)d^(p{x), (E.l) 


This implies that G{x, y) can only be a function of quanti¬ 
ties that transform as scalar bilinears. An obvious example 


is the geodesic distance s(x, y) between the two points x 
and y, 


"^{v) 

, \ f ,~ds{T) 2 - dx^^{ t) dx''{f) 2 

six,y) = J dr-^ , ds =9^.^ - ^dr , 

f(x) 

Ss{x,y) = - (E.2) 


with x^{t) on a geodesic parametrized by an invariant pa¬ 
rameter f. The propagator in flat space (EZD has indeed 
the form 


G{x,y) 


1 


(E.3) 


For a geometry with curvature there are other possible 
quantities transforming as scalar bilinears. For example, 
we may replace ds by ds where in the definition of is 
replaced by or in the definition of ds. For the spe¬ 
cial case of a maximally symmetric space all geometric ten¬ 
sors are constructed from the metric, ^ R =const, 
etc.. Then s is the only independent bilinear function of 
X and y, in the sense that all other covariant bilinears are 
functions of s. This applies to de Sitter space. A Green’s 
function that only depends on geometrical quantities in 
de Sitter space (and not on boundary conditions or other 
fields) can therefore only depend on s, G(x, y) = G(s). 


2. Global S'0(1,4)-symmetry of de Sitter space 

If the background geometry admits isometries a unique 
scaling correlation must be invariant under the correspond¬ 
ing symmetry transformations. This is supposed to be re¬ 
alized as a consequence of the time evolution even if the 
initial correlation function does not exhibit this symmetry. 
Invariance under symmetry transformations constitutes an 
important restriction for the allowed function space for the 
Green’s function. The geodesic distance is invariant under 
such symmetry transformations, e.g. the Poincare trans¬ 
formations for flat space or S'0(l,4) transformations for 
de Sitter space. For G{x, y) = G{s) the invariance of the 
propagator is guaranteed. 

For de Sitter space the four generators of the symmetry 
group SO{l, 4) beyond the three dimensional rotations and 
translations are 

Mo = gdr,+x^dj, (E.4) 

A'h = Ho [2xigdr, + 2x^x^dj + [if - |^|^)5i] . 

A scalar bilinear G(a;, y) transforms under the correspond¬ 
ing infinitesimal transformation as 

5G{x, y) = {a[Mo + Mo) + PiiM^ + Mi)}G{x, y), (E.5) 

with M obtained from M by y —>■ 77 ' Xi —>■ j/i . For functions 
G(r,r],g') depending only on r one has 

Mo + Mq = rjdn + T]'drf'+ r6r, 

Mi Mi = Ho[xi yi){Mo Mo) Hoixi — yi)Mr, 

2 /2 

Mr = gdn — g'dn' + - -^9^.. (E. 6 ) 
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An 4 ) invariant Green’s function must therefore obey 

{Mq + Mq)G = 0 , M^G = 0 . These relations hold for the 
combination 


we find that the first term in Gg in eq. ( 1821 ) is invariant 
provided that a does not depend on k. For the other pieces 
we compute (/ 3 + = (/3 + *7)/2) 


— {rj — rj'Y 

m' 


(E. 7 ) Mo [^+(fc)u;+(? 7 )u;+(??')]=(»?'). (E. 15 ) 


e.g. (Mo + Mo)z = 0 ,MrZ = 0 . The invariant 0 is related 
to the geodesic distance in de Sitter space [ 1 ^, @ 


z = 2 {l-P) , P = 


cos(Mos) for > 0 

cosh(jiJos) for < 0 ' 


The propagator in de Sitter space is therefore only a func¬ 
tion of z, G{r,r],ri') = G{z). 


3. De-Sitter symmetry in momentum space 

The action of the symmetry generators on scalar bilin- 
ears in Fourier space can be inferred by inserting eqs. (SHI), 
SH) into eq. (lE^SI) . 


(Mo + Mo)Gik, r,, r,') = M^G{k, r,, rj'), (E. 9 ) 


and 


(M, + M,) J P^^^-y^Gik, 77,77') 

k 

= HoJ {{x, + 7 /,)MoG(fc, 77,77') 

k 

+ 7 ^M,G(fc, 77 , 77 '), (E.IO) 


with 


d , d , 9 „ 

Mq = ^ Try - - 3> 

arj arj' ok 


d 


, d 


d 




dk 

An invariant correlation function has to obey 
MoG{k, T], 77') = 0 , MrG{k, 77,77') = 0 . 


(E. 12 ) 


We observe that the symmetry transformations mix differ¬ 
ent fc-modes due to the appearance of d/dk in the genera¬ 
tors. 

We may now discuss the symmetry transformations 
of the different pieces of the general propagator in eqs. 
(|M 1 ) . (IM)) with given by eq. m . One finds that Ga is 
invariant 

Mo [wfe {v)wt iv')] = 0 , Mr [w^ (77)77;+ (77')] = 0. (E. 13 ) 

(The invariance of (7)^^ (77') follows by complex conju¬ 
gation.) With 


Mo [a{k)w^ (77)77;+ (77')] = iv') (E. 14 ) 


such that symmetry requires a fc-independent constant /?+. 
With 


M,. [77;+(77)77;+(77')] = 0 (E. 16 ) 

we conclude that the combination (??)'w^ (77') is also in¬ 
variant under the symmetry of de Sitter space. By virtue 
of complex conjugation the same holds for 777^(77)77;^ (77'). 
We conclude that for /c-independent a, (3 and 7 the gen¬ 
eral symmetric part Gg in eq. ( 1821 ) is invariant under the 
de Sitter symmetry. In contrast to Poincare symmetry for 
flat space, de Sitter symmetry does not fix the Green’s 
function uniquely. One remains with a family of invariant 
propagators, parametrized by the three constants a, f 3 and 
7 - 

4. Scale invariance 

The generator Mq corresponds to global scale transfor¬ 
mations of coordinates. A scale invariant correlation func¬ 
tion obeys 


(Mo -h Mo)G = 0 . (E. 17 ) 

In momentum space this simply means that kz’G can only 
depend on the dimensionless combinations krj and krj'. 
From eqs. (|E. 14 p . (IE. 15 l) it follows directly that scale sym¬ 
metry alone implies that a, /3 and 7 are independent of k. 
Indeed, the combinations k^w^{ri)w^{r]') already depend 
only on kr] and kr]' and are therefore scale invariant. Any 
fc-independence of a, / 3 , 7 would have to come along with 
a dependence on 77 or 77' which is contradicting their role 
as integration constants. 

Since for fc-independent a, / 3 , 7 the correlation func¬ 
tions are also independent under the full 5 ' 0 (I, 4 ) de Sit¬ 
ter symmetry, one concludes that scale symmetry, together 
with the symmetry properties and the propagator equa¬ 
tion for a massless scalar, implies full S' 0 (l, 4 ) invariance 
of the correlation function. (One should notice that this 
is not a pure symmetry property. There are many func¬ 
tions that depend only on kr] and kr]' and they all would 
lead to scale invariant k^G. Only very particular ones are 
S' 0 (l, 4 ) symmetric, however.) The impact of scale invari¬ 
ance extends to the propagator of a massive scalar field in 
de Sitter space. Scale symmetry implies that coefficients of 
the combinations A:^77;+ (77)77;+(77') must be fc-independent, 
if 77;+ = 7 ;+/a are the two independent solutions of the 
homogeneous equation DiP^ = 0, normalized according to 
ea. (II 34 l) with 7;+ = v^{y)/V^- Indeed, for de Sitter space 
one has y = —krj such that k^w^ {ri)w^ (rj') is only a func¬ 
tion of krj and ki]'. No further fc-dependence of coefficients 
is allowed by coordinate scale symmetry. 

We may recall here the different status of coordinate 
scale symmetry and dilatation symmetry. Coordinate scale 
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symmetry is part of the general coordinate transformations 
with = ex^. The effective action is therefore invariant. 
This holds in the presence of arbitrary couplings with di¬ 
mension of mass, as, for example, a scalar mass term. Co¬ 
ordinate scale symmetry can be broken, however, by the 
geometry and field values of a given cosmological solution. 
On the other hand, dilatation transformations are rescal¬ 
ings of fields according to their (canonical or anomalous) 
dimension. It is typically violated by couplings with di¬ 
mension of mass since they single out a particular scale. 
Dilatation symmetry of the effective action restricts its 
form far beyond general coordinate invariance. The con¬ 
nection between the two symmetry transformations is that 
the field equations of a dilatation invariant effective action 
may have solutions that preserve dilatation symmetry (no 
spontaneous dilatation symmetry breaking). For such so¬ 
lutions coordinate scale symmetry is realized. 


5. Propagator in de Sitter space 

According to our conjecture the scaling correlation for a 
massless scalar field in de Sitter space has in momentum 
space (for fc ^ 0) the unique form ea. d^ . The lack of 
regularization leads for fc —^ 0 to an infrared behavior that 
does not admit a Fourier transform. In position space no 
propagator consistent with our conjecture and the symme¬ 
tries of de Sitter space exists. On the other hand, one can 
find explicit solutions of the propagator equation in posi¬ 
tion space which are consistent with the symmetries. Since 
the issue of the correct choice of the propagator is debated 
(20l - [^ Ibbl - f^ we sketch here a direct computation. This 
will reveal the connection to the general form of the propa¬ 
gator in momentum space. The problem of the absence of 
a propagator in position space for our conjectured unique 
scaling correlation disappears if one regularizes the inverse 
propagator by adding a small positive mass term. 

In position space the propagator can only depend on the 
invariant z, given by eq. (lEJl) . or the associated invariant 


z rf' rj''^ — 

2 2?7?7' 


(E.18) 


For functions G{P) depending only on P one has 


DG=^bpG, 


(E.19) 


with 


bp = (p2 _ i)dl + dPdp. (E.20) 

If we add a mass term for the scalar field, m? = rn?/H q^ 
the defining equation for the Green's function reads 

{bp + mbG{P) = -iH^S{r)S{r] - r]'). (E.21) 

Let us first consider for r ^ 0 or rj ^ r]' the homogeneous 
equation 

[{P"^-l)dl + 4Pdp + m'^]G{P) =0. (E.22) 

This equation is symmetric under the replacement P —>■ 
—P. It is convenient to employ 

G=(l-p2)-5P, (E.23) 


such that the homogeneous equation becomes 

I (1 - Pbdl - 2Pdp + 2 - P(P) = 0. 

(E.24) 

This is the defining equation for the associated Legendre 
functions P^ and Q},, with 

p{p + l) =2-ffb. (E.25) 


(The index iz of the associated Legendre functions should 
not be confounded with the symbol p used in the preceding 
section for H/H^.) 

The general solution of ea. (IE.22|) is therefore given by 

G = (1 - P2)-^ {B^PliP) + B^QliP)} . (E.26) 

For a massless scalar field {fh^ = 0) one has p = 1 and the 
general solution of the homogeneous equation reads 


G = Bi+B2 


P 

1 - p2 


-\n 

4 


/( 1+^\1 


(E.27) 


corresponding to the proposal in ref. [l^, see also 
This solution has two poles at P = 1 and P = —l,orz = 0 
and z = 4. Comparing the behavior near P = 1 with the 
free propagator yields P 2 = Hq/{A'k"^), or 


r = ^ 

47r2 


1 1 , 

- -I- - In 
z 4 


(4-z)^ 


4- z 


Pi. (E.28) 


For Pi = 0 this is identical with G in eq. (IA.22I) . We 
observe G(z = 2) = Pi. 

We next turn to the inhomogeneous equation 

+ (2 -z)9^| G= ^^^^d(r)d(ry- 77 '). 

(E.29) 

For G the pole at z = 0 generates the r.h.s. of ea. (IE.29|) . 
fixing the normalization (IE.281) . The second pole at z = 4 
or = {r] + produces a similar inhomogeneous term 
at T] + r]' = 0. For the allowed range of 77 < 0, ry' < 0 for de 
Sitter space this inhomogeneous term is not encountered. 
We conclude that G obeys the propagator equation. In the 
following we set Pi = 0 such that eqs. (IE.28|) and (IA.22|) 
coincide. The form of G in momentum space is given by 
eq. (IA.21I) . The propagator G is not compatible with our 
conjecture for the scaling correlation. 

We have already discussed the power spectrum of pri¬ 
mordial fluctuations for the propagator G in sect. IIVI It 
differs from the one of the Bunch-Davies vacuum by a sup¬ 
pression of the overall amplitude. There seems to be no 
obvious observational possibility to distinguish between G 
and the de Sitter propagator (1^ . The question which one 
of the propagators is “correct” is finally an issue of the time 
evolution of correlation functions as discussed in sect. IS 
We observe that the form of G has its own diseases. The 
equal time correlation function turns negative in the re¬ 
gion < 477 ^ (z < 4) and has a pole for = 477 ^. On 
the other hand the absence of a Fourier transform for the 
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propagator dsni) disappears in the presence of a small scalar 
mass term. 

Indeed, the non-existence of the correlation function in 
position space is directly linked to the lack of invertibility 
of the inverse propagator For a massive scalar field, 

> 0 , the inverse propagator has no longer a vanishing 
eigenvalue and the propagator exists both in position and 
momentum space. Adding an arbitrarily small mass term 
regulates the propagator. We have already discussed 
this propagator in momentum space in appendix D. For 0 < 
rh^ <C 1 the behavior of the equal time Green’s function 
for small k is given by 


lim G(fc, f]) 
fe—>-0 


2fc3 \aHj 


(E.30) 


for G{P), and the normalization corresponds to a leading 
behavior G{P —?> 1) = iF^/( 87 r^(l — P)). A cancellation of 
the pole at P = — 1 is not possible for = 0 since one 
of the basis functions for the most general homogeneous 
solutions is simply a constant. 

6. Special case iri^ = 2H^ 

An interesting special case for a massive propagator is 
= 2. In this case the general solution of the homoge¬ 
neous equation ljE.2211 reads 

= + (E.31) 

and the inhomogeneous equation (jE.2I|) is obeyed for o_ = 

0 , 


The Fourier transform G(r, rj) is well defined. Since the 
mass term respects the symmetries of de Sitter space 
G(r, 77 , 77 ') is only a function of z. It is a hypergeometric 
function, corresponding to the linear combination (IE.261) 
with u = 1 — rh^j3. After imposing the normalization 
via the inhomogeneous term and the reality condition we 
remain with a family of de Sitter invariant propagators, 
labeled by a complex constant. They are associated with 
the “a-vacua” [lo, [ 2 l|. 

For the scaling correlation the coefficients Bi and B 2 
should be chosen such that no pole at P = — 1 occurs 


87r2(P— I) 

477^0(77)0(77') — (77 — 77')^] ° ^ ^ 

Indeed, the time evolution eouation d 17211 simplifies to 

{dn^ + k'^)v = 0, (E.33) 

which is the same as for flat space. 
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